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SIITQPSIS 

OPTIMAL ATTITUDE STABILIZiTIOU AED OOMTROL OF SPAGECRAFT BY 

SOLAR RilDIATION PRESSURE 

A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 
DOCTOR OF PHILOSOPHY 

by- 

Rapur Y enkat achalam 
to the 

Department of Mechanical Engineering 
Indian Institute of Technology, Kanpur 

May, 1980 

Attitude control of high altitude spacecraft by. means of 
solar radiation pressure is investigated. The objective is to 
develop solar pressure attitude control sy^ems that are 
mechanically simple to implement and yet offer desirable control : 
characteristics. A solar controller consisting of two highly 
reflective lightweight control surfaces rotatable about body- 
fixed axes is proposed for a variety of spacecraft missions. 

Such a configuration represents the minimum hardware 
implementation as at least two control surfaces are required to 
achieve the solar torque reversals. Spacecraft missions 
demanding earth-pointing attitude, inert! ally -fixed orientation 
and large-angle attitude manexivers in presence of gravity- ; 

gradient torques are considered. Control laws directly governing 
the rotations of the two control surfaces are synthesized to 
meet the mission requirements and their validity is established 
for all times of the year. 
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The problem of stabilizing the pitch attitude of an 
unsymmetrical earth-pointing spacecraft is formulated first, 
iilthough the pitch dynamics is relatively simple, the solar 
radiation torque generated by the controller turns out to be 
a rather complicated transcendental function of the pitch 
attitude, control surface rotations and the time. To make 
the control synthesis mathematically tractable, linearization 
of the governing nonlinear equation of motion for small amplitude 
pitch motion and control surface excursions is undertaken, A 
piecewise linear feedback control law with constant gains is 
proposed for the differential rotation of the control surfaces. 
The system stability as a function of the attitude error and 
rate gains is established using T^loquet^ theoiy. The approach is 
next extended to examine the potential of the controller in 
achieving three-axis Stabilization of unsymmetrical earth- 
pointing spacecraft. Optimal Control theory is employed to 
synthesize a feedback control law with periodic time-varying 
gains for the differential control surface rotation leading ■ 
to an asymptotically stable system. Subopt imal control 
strategies involving gain scheduling a few times per year are 
suggested which lead to considerable reduction of the associated 
computational requirements. 

Applications of the two-surface solar controller for 
missions demanding stabilization along arbitrary inert ially-fixed 
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attitudes are considered next. The inertial pointing problems 
are characterized by the need to provide nominal controls 
countering the time-varying gravity gradient disturbances. 

Por pitch stabilization, both the nominal and feedback 
control laws governing the differential control surface 
rotation are synthesized analytically. The stability properties 
of the system are shown, through certain transformations, to be 
governed by the classical Mathieu equation. The investigation 
is then extended to the problem of imparting inertially -fixed 
orientations to the symmetry, axis of an axisymmetric spacecraft, 
Por this more complex problem, a numerical approach is followed. 
Nominal control laws for the individual control surface rotations 
are developed and the governing equations of motion linearized 
in the neighbourhood of the desired inertial attitude. 

Minimization of a quadratic performance index then results in an 
optimal feedback control law for the differential control surface 
rotation assuring asymptotic stability for the system. Suboptiraal 
control policies promising computational reductions during system 
implementation are explored. 

Finally, the feasibility of achieving large-angle attitude 
maneuvers of the spacecraft symmetiy axis is investigated. The 
attitude transitions being unrestricted, the full nonlinear 
equations of motion in the gravity-gradient field are treated 
for the ^nthesis. The problem is formulated as a minimum time 



transfer problem with terminal state constraints. A gradient 
procedure is developed to determine the optimal open-loop 
strategy governing the differential control surface rotation 
for the specified maneuver. The transfer strategy is then 
coupled, with the appropriate nominal and feedback control 
policies required for stabilization along the inertially- 
fixed terminal orientation. 

The control system performance is analyzed through 
simulations of the complete nonlinear equations of motion in 
both eairfch-pointing as well as inertially-fixed modes of 
stabilization. The loss of solar control during the earth’s 
shadow passage by the spacecraft is also modelled. 

Disturbances representing micrometeorite impacts and initial 
attitude misalignments are introduced. The results indicate 
a moderately sized controller to be capable of achieving 
control times of the order of a fraction of an orbit. 

The thesis establishes the effectiveness of a 
mechanically simple two-surface solar pressure control system 
in stabilizing the spacecraft along earth-pointing as well 
as inertially-fixed attitudes. The analysis demonstrates 
the capability of the controller in stabilizing the spacecraft 
even along gravitationally unstable equilibria. Its ability 
to achieve stabilization along arbitrary inertially-fixed 
orientations opens new vistas for the concept of solar 
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pressure attitude control. The feasibility of accomplishing 
large-angle attitude maneuvers adds versatility to the concept 
enabling the spacecraft to undertake diverse missions. The 
approach is semipassive -in 'Character as no mass expulsion 
schemes are involved and only a small power consumption is 
envisaged to rotate the control surfaces. This promises a 
long operational life-span for the spacecraft. 
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1. IMDRODUO!I?IOir 
1.1 Preliminary Remarks 

Phe general motion of a space vehicle consists of the 
translation of its center of mass and the rotational motion about 
the center of mass. The former describes the orbital motion of 
the vehicle and the latter is referred to as its attitude motion. 

The attitude dynamics of a satellite is subject to a 
variety of environmental influences. These are ; the gravity- 
gradient effect, the earth's magnetic field, aerodynamic forces, 
radiation pressure due to the earth's radiation and albedo at 
low altitudes ( < SOO km), and direct solar radiation pressure 
at high altitudes. Attitude motion may also be excited due to 
encounters with micrometeorites. On the other hand, successful 
operation of a vast majority of space missions demands that the 
spacecraft maintain a preferred orientation in space. Por 
example, communications, weather, military, earth-resources 
applications, etc., require the spacecraft to continually point 
towards the earth. Missions involving scientific and astronomical 
observations may demand the satellite to maintain an inertially- 
fixed pointing. The capability to accomplish attitude maneuvers 
between any two specified spatial orientations may also be 
desired in certain applications. An effective attitude control 
system is hence necessary to achieve the mission objectives. 
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1. IMDROJDUOTIOJr 
1.1 Preliminary Remarks 

The general motion of a space vehicle consists of the 
translation of its center of mass and the rotational motion about 
the center of mass. The former describes the orbital motion of 
the vehicle and the latter is referred to as its attitude motion. 

The attitude dynamics of a satellite is subject to a 
variety of environmental influences. These are ; the gravity- 
gradient effect, the earth's magnetic field, aerodynamic forces, 
radiation pressure due to the earth's radiation and albedo at 
low altitudes ( < 800 km), and direct solar radiation pressure 
at high altitudes. Attitude motion may also be excited due to 
encounters with micrometeorites. On the other hand, successful 
operation of a vast majority of space missions demands that the 
spacecraft maintain a preferred orientation in space. Por 
example, communications, weather, military, earth-resources 
applications, etc., require the spacecraft to continually point 
towards the earth. Missions involving scientific and astronomical 
observations mqy demand the satellite to maintain an inertially- 
fixed pointing. The capability to accomplish attitude maneuvers 
between any two specified spatial orientations may also be 
desired in certain applications. An effective attitude control 
system is hence necessary to achieve the mission objectives. 
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Several methods are available for attittide control which 
involve the use of active control elements, spin-stabilization 
technique and the utilization of environmental forces. Active 
stabilization procedures involve mass expulsion schemes and/or 
components requiring a large amount of energy which is an 
expensive commodity aboard a spacecraft. The advantage of the 
technique lies in the ability to achieve any orientation with a 
very high degree of accuracy. However, the large energy 
consumption leads to a reduced operational life-time for the 
satellite. Spin-stabilization relies on the inherent tendency 
of a spinning body to maintain a fixed attitude in torque-free 
space. The method does not involve any energy consumption. Its 
pointing accurac}'', however, is limited as it deteriorates due 
to the effect of external disturbances. Combinations of the 
active methods and spin-stabilization are commonly employed in 
the form of momentum bias or dual-spin configurations. 

Methods of attitude control based on the utilization of 
the environmental forces are either completely passive or 
involve very small energy consumption. The gravity-gradient 
stabilization technique is based on designing the satellite 
mass distribution such that the earth-pointing orientation of 
the satellite represents a stable configuration in the gravity- 
gradient field. Although the procedure involves no energy 
consumption, low pointing accuracy and very slow rates of 
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attitude correction are achieved. The exploitation of the 
aerodynamic forces, the earth's magnetic field and solar 
radiation pressure involve the introduction of suitable control 
elements which may be operated with very little electrical power 
to generate attitude control torques. Por example, movable 
aerodynamic surfaces m^ be employed at low altitudes to 
intercept the relative free-molecular flow. Onboard 
electromagnets are used to generate control torques through 
interaction with the earth’s magnetic field over a range of 
orbital altitudes. For high altitude spacecraft, control 
surfaces are employed to intercept the incident solar radiation 
to develop control moments about the satellite center of mass. 
These methods promise high pointing accuracies with an 
increased operational life-span for the spacecraft. 

This thesis deals with the utilization of solar radiation 
pressure for the attitude control of high altitude spacecraft. 
The development and analysis of solar pressure attitude control 
sy stars for a variety of spacecraft missions forms the main 
objective of the work. 

1,2 literature Review 

Most of the early researches in the utilization of solar 
radiation pressure for attitude control relate to the problem 
of maintaining a sun-pointing spacecraft attitude. Sohn (1959) 
considered solar pressure attitude control of an interplanetary 



4 


spacecraft using a weatheirvane type tail surface. This enables 
the spacecraft to maintain a nominally fixed orientation relative 
to the satellite-sun line. The author also mentioned the idea of 
including rate sensors for solar damping of undesirable 

2 

oscillations about the nominal orientation, Prye and Stearns 
( 1959 ) presented a general survey concerning the utilization of 
gravity-gradient, aerodynamic and solar radiation torques for 
attitude stabilization and also provided a brief description of 
sensors and control devices. The authors proposed a trailing 
cone system to intercept the direct solar radiation to derive 
control torques tending to align the cone axis with the sun-line. 
The concept of introducing small servo driven solar control 
surfaces, called ’Solarons', to achieve attitude damping was 
also mentioned. ITewton (i960) suggested the possibility of 
generating corrective radiation torques by designing spherical 
satellites with one-half of the surface having reflective 
properties and the other half absorptive properties. ¥ith the 
reflecting half nominally facing the sun, a stabilizing 
radiation torque is achieved. Hibbard^ (1961) proposed a 
reflector-collector configuration for providing focussed 
radiation pressure with a view to achieve onboard weight 
savings. Merrick et al.^’ (1966, 68) considered the damping 

of attitude oscillations by providing relative motion between 
a pair of connected bodies. Acord and SFicklas' (1964) evolved 
a passive stabilizer capable of providing attitude damping 
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through a thermo-mechaaical phase lag arrangement. 

The problem of aligning the spin-axis of a spinning 

8 

spacecraft towards the sun was investigated by Ule (1963). 

The author suggested the use of windmill type mirror-arrays 

and corner mirror-arrays for producing restoring precession 

torques. A mechanism for spin rate control involving a pair of 

flat surfaces with one side reflective and the other side 

absorbent was proposed. The angle of attack of the spin rate 

control surfaces could be governed either by remote control or 

by centrifugal forces. A passive solar damping approach where 

a structural element provides a phase shift between the solar 

restoring torque and attitude motion was suggested by Carrel 

9 

and Limburg (1965). Purther development of various mechanisms 

10 1 *1 

along similar lines was done by Peterson (1966), Colombo 

1 2 

(1966) and Harrington (1966). The use of a grated solar sail 

for attitude control of sun-pointing spinning spacecraft was 

1 1 4 - 

studied by Palcovitz (1966). Por the same problem, Crocker 

(1970) proposed two body mounted differentially coated solar 

surfaces for precession control and a pair of spring mounted 

15 

solar paddles for spin rate control. Pande (1976) considered 
semipassive spin-axis orientation control by means of two 
rotatable control surfaces. The concept enables the spin-axis 
to attain any arbitrary pointing in space. 

The above mentioned studies represent rather simplified 
situations as a torque-free environment is assumed throughout. 
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Solar pressure attitude control of satellites under more 

i S 

realistic dynamic conditions was examined by Mallach (1966). 
The author suggested the possibility of achieving solar 
damping of the librations of satellites in the gravity- 
gradient field and presented a simplified analysis using 
averaged radiation torques. Modi and Planagan (1971) 
examined the planar attitude control of a gravity -gradient 
system in an ecliptic orbit using the radiation pressure to 
provide the damping torque. Modi and Tschann * ^ (1971 ,73) 
extended the analysis by the introduction of a displacement 
and velocity sensitive controller. Practical feasibility of 
the system was considered through an arrangement involving an 
unfurl able membrane whose effective area could be judiciously 
governed. Modi and Kumar ’ (1971,74) studied planar as 

well as coupled libration control through the use of solar 
surfaces which may be translated relative to the spacecraft 
body to govern the solar control torques, 

A more practical approach, employing highly reflective 

control surfaces rotatable about body-fixed axes, was proposed 

22 

by Modi and Pande (1973) for three-axis control of a dual- 
spin satellite, A large number of control surfaces were found 
necessary to achieve independent control of the roll, yaw and 
pitch torques. Unf ortunately , the solar torque components turn 
out to be extremely complicated, coupled, transcendental 
functions of the control surface rotations and the orbital 
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angle. Therefore, a coupled system of transcendental equations 
must be solved continuously to determine the required control 
surface rotations. This limits the applicability of the method 
to spacecraft having a sophisticated onboard computational, 
capability. Modi and Pande (1974) also developed a simpler 
model involving four independently rotatable control surfaces 
for three-axis control of axisymmetric spacecraft. Bang-bang 
control laws directly specifying the control surface rotations 
as functions of the attitude errors and rates were suggested. 
This not only led to computational simplifications but also 
improved the rates of attitude correction for large errors. 
However, an "undesirable chatter behaviour was observed in the 
small epror regime. The analyses investigated the applicability 
of the technique in nonecliptic orbits of arbitrary inclination 
and suggested modifications to the control laws to reduce the 
eccentricity induced steady state errors, Pande et alt ^ (1974) 
investigated the minimum time pitch control problem for a 
gravity-orient ed satellite. The authors also recognized the 
possibility of using solar pressure control in conjunction with 

pir p^ 

magnetic (1974) and aerodynamic (1974) modes of control. 

There have been two experiments to date to examine the 
practical performance of solar pressure attitude control 

systems. The first relates to the flight of Mariner 17 

27 2R 

spacecraft ^ (1966, 69) which was equipped with four solar 
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vanes attached to the tips of the solar panels. The vanes 
were designed to passively stabilize the spacecraft so as to 
align its roll axis along the snn-line. A small amonnt of 

7 

damping was also provided hy the passive solar vane actuator. 
Unfortunat ely , one of the vanes proved to he inoperative 
during the early portion of the mission. Subsequently, a power 
transient aboard the spacecraft caused reactivation of the 
solar vane. The solar pressure experiment, conducted on 
depletion of the attitude control gas, showed the solar vanes 
to be capable of maintaining the spacecraft roll axis within 
1'^ of the sun-line in conjunction with nutation damping 
provided by gyros. A recent experimait ^ (1979) conducted by 
the European space Agency established the effectiveness of the 
method for geostationaiy communications satellites. The solar 
panels of OTS-2 spacecraft were rotated differentially to 
generate control torques about the roll and yaw axes. A ground 
commanded discontinuous mode of control was used for the solar 
panel rotation during the experiment. For a period of six days, 
all the spacecraft thrusters were completely disabled and its 
attitude was controlled successfully entirely by solar radiation 
pressure, 

1,5 Purpose and Scope of the Investigation 

From the foregoing, it is apparent that most practical 
solar controller conf igurations proposed for earth-pointing 
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spacecraft involve a nnmPer of rotatable control surfaces. 
Clearly, it -woull be of interest to simplify the hardware 
implementation through the minimization of the number of control 
surfaces. Moreover, the continuous feedback control laws 
developed for the solar torque components require extensive 
onboard computational capability to determine the commanded 
control surface rotations. On the other hand, control 
strategies directly specifying feedback relations for the 
control surface rotations are limited to a discontinuous mode 
of operation leading to an unsatisfactory perfomance in the 
small error region. Hence, there is a need to develop continuous 
control laws that directly govern the control surface rotations 
as functions of the attitude errors and rates. 

The literature review also indicates that utilization of 
solar pressure for attitude control has been limited to either 
sun-oriented or earth-pointing spacecraft. Its application for 
attaining space-oriented attitudes, unfortunately, has remained 
unexplored. On the other hand, solar pressure stabilization of 
satellites along an inert ially-fixed attitude would represent 
an inexpensive technique for long-life applications such as 
astronomical observations, furthermore, it would be of interest 
to develop solar controllers capable of achieving arbitrary 
large-angle attitude maneuvers in space. Such a system would 
add versatility to the spacecraft and enable it to undertake 
diverse missions through an entirely semipassive approach. 



The thesis makes aa attempt at developing solar pressure 
control systems meeting the above mentioned objectives. 
Controller configurations consisting of two highly reflective 
lightweight control surfaces rotatable about body-fixed axes 
are proposed for the different pointing objeotives. Such 
configurations represent the minimum hardware implement ation 
as at least two control surfaces are required to achieve the 
control torque reversals. 

The problem of stabilizing the pitch attitude of an 
unsymmetrical earth-pointing spacecraft is analyzed first. 

A piecewise linear feedback control law with constant gains is 
proposed for the differential rotation of the control surfaces. 
The system stability and performance are analyzed as functions 
of the attitude error and rate gains. 

Fext, the approach is extended to develop a solar control 
system for three-axis stabilization of unsymmetrical earth- 
pointing spacecraft. Optimal control theory is employed to 
synthesize a feedback control law for the differential control 
surface rotation leading to an asymptotically stable system. 
Suboptimal control strategies are suggested which promise 
considerable reduction of the computational requirements during 
the system implementation. 

This is followed by an investigation of missions demanding 
an inertially-fixed attitude. The inertial pointing problems 
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are oharact erized by the need to provide nominal controls 
countering the time-varying gravity-gradient disturbances. 

Por pitch stabilization, both the nominal and feedback control 
laws governing the differential control surface rotation are 
synthesized analytically. The influence of system gains and 
parameters on its stability and response characteristics is 
established. 

The study is then extended to the problem of imparting 
inerti ally-fixed orientations to the symmetiy axis of an 
axisymmetric spacecraft. Por this more complex problem, a 
numerical approach is followed. iTominal control laws for the 
individual control surface rotations are developed and linear 
optimal control theory is applied for the control synthesis in 
the neighboirchood of the desired inertial orientation. Again, 
the possibility of employing subopt imal control policies with 
reduced computational requirements is explored, 

Pinally, the feasibility of achieving large-angle attitude 
maneuvers of the spacecraft symmetry axis is investigated. 

The problem is formulated as a minimum time transfer problem 
with terminal state constraints. A gradient procedure is 
developed to determine the optimal open-loop strategy for the 
differential control surface rotation. The transfer strategy 
is then coupled with the appropriate nominal and feedback 
control policies required for stabilization along the inertially 
fixed terminal orientation. 

5'igure 1.1 schematically represents the plan of study. 



Optimal attitude stabilization 
and control of spacecraft by 



Fig. II Plan of study 
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2. PITCH iTTITDDE STABILIZATION OP EAREH-POIUTING SPAGEGRAPT 

This chapter presents the development of a solar controller 
for controlling the pitch attitude of an unsjrrametrical earth- 
pointing satellite. A controller configuration involving two 
highly reflective rotatable control surfaces is proposed. The 
governing equation of motion is found to be nonlinear and 
nonautonomous in character, rendering it difficult to arrive 
at a suitable control strategy for the control surface rotations. 

A simple feedback control law, directly specifying the differential 
rotation of the control surfaces, is developed following 
linearization of the pitch equation. The stability of the control 
system is analyzed using Ploquet theory and the validity of the 
control law is established for all times of the year. The system 
performance is studied through simulation of the full nonlinear 
equations of motion. 

2,1 Problem Pormulation 

Pigure 2.1 shows the geometiy of the orbital and attitude 
motion of an unsymmetrical satellite with its center of mass S 
moving in a circular orbit about the eaih;h's center 0. The 
inertial frame of ref erence,X,Y,2,is selected such that the 
Y-axis points towards the ascending node N, and YZ defines the 
plane of the orbit. The unit vector e indicates the orientation 
of the earth-sun line, x^, y^, z^ represent the rotating 
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ortital coordinate system with y^ along the local veitical, 

Zq tangent to the orbit in its plane and along the orbit 
normal. The spacecraft executes pitch attitude motion \ in 
the plane of the orbit. 

The solar pressure controller consists of two identical, 
highly reflective, lightweight surfaces and P 2 "which are 
permitted rotations *5^ and 52» ^respectively , from the y-axis. 
The center of pressure of each surface lies on the satellite 
y-axis (may be anywhere in yz-plane). This, along with the 
fact that the radiation force on a highly reflective surface 
is directed along the surface normal, ensures that only a 
pitch moment is pixjduced by the controller. 

The equation of pitch motion may be written as 

(d/dt) tl ( X + e) J = (2.1) 

where M_ and represent the gravitational and the solar 

& s 

torques about the pitch axis. 

The gravitational torque component M is well known to 

r, 0 

w JO 
be , 

= -3(1^ - I,Jsin X cos X (2.2) 

Pigure 2.2 shows the components of the radiation force 
acting on the control plate P^^ under specular reflection 
conditions. The j resultant solar pressure force acting on 
the plate is gk pa by 



Normal 



Fig. 2-2 Radiation forces acting on control 
plate Pj 
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= -pAjcos I C(1-T)e + P (2.3) 

Expressing 'a^ in terms of e and n^, and ignoring terms of 
order (1-T-p)/2p compared to iniity, which is justifiable for 
surfaces of high reflectivity, it simplifies to 

= -2p pA I cos cos (2.4) 

The angle of incidence Cj_ is obtained from 

cos = e.n^ (2.5 ) 

The unit vector along the sun-line may be expressed as 


e = 1 + By 3 + k (2.6) 

Where 

e = - sin <P sin i (2.7a) 

e = cos cos(9 +X ) + sin <l> cos i sin(9+X) (2.7b) 

V 

Sg = -cos 4» sin(0 +x ) + sin<J>cos i cos(9+x) (2.7c) 

The unit normal to plate is given by 

= sin J - cos S' (2,8) 


Substitution of Equations (2.6) and (2.8) in Eqiation (2.5) 
gives 

2 2 4 - 

cos = (1-sin sin i)"^ sin(e+ 5 + X + 6 ^) (2.9) 

where 

C (4>) = -arctan(tan <j> cos i) 


( 2 . 10 ) 



The solar pressure moment acting on the spacecraft then becomes 

Mg = c J Z - el X fg (2.11) 

= 2 p pA e ( j cos j cos cos 6 ^ - 

j cos ^2 I 5 2 ^ ( 2 . 12 ) 

Substituting for M and M from Equations (2.2) and (2.12) in 

g s 

Equation (2.1), and changing the independent variable from t 
to 6 through the relation e = flt, the governing pitch equation 
takes the form 

X” +32 sin X cos ^ (2.13) 

where 

Mg = CaClsin(0+5+X+5^) I sin(e+ 5 +X+ 6 ^) cos 5^ 

- j sin(9+c+X+ 62^1 sin( 0 +?+x+ 62 ^°°® ^2^ (2.14) 

Here, the function a (^ ) is given by 

a ((f») = 1 - sin^<|> sin^i (2.15) 

and the solar parameter , C, characterizing the controller size 
is defined as 

0 = 2 p pAe /I (2.16) 

As varies by 2 ir radians in a year, both ©(4' ) and C (<f>) 
represent extremely slow variables compared to the orbital 
angle e. Therefore, they remain practically constant over a 
number of the satellite orbits. In the analysis to follow, 
the solar aspect angle 4> as well as other slow variables such 
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as cy(4>), ?(<^)} etc., are treated as system parameters which 
acquire different values at different times of the year. 

2.2 Control Synthesis 

When the two control surfaces are parallel (6^ = ^2 = 6g) 
the stable equilibrium position of the system corresponds to 
the alignmoat of the minimum moment of inertia axis along the 
local vertical, i.e. , x = 0, ir with K > 0. Consequently , a 
control law of the form 

62,^ ,e ) = -11 X’ - VX ( 2 . 17 ) 

where P , v represait constant positive gain factors, leads 
to asymptotic stability of the pitch motion* However, its 
implementation would require continuous solution of the 
transcendental Equation (2.14) to determine the 5^, 62 
required to generate the commanded torque, iilso , as this 
would mean solving one equation for the two unknowns, a 
suitable constraint between them must be specified. 

In order to bypass the computational difficulty, a 
simplification of the torque expression is presently 
undertaken. It is assumed that the control surface excursions 
from their ’off’ positions are small. Replacing 6^ by 
(&Q + ^2 by (5^ + retaining only the linear 

terms in 5^ and 62’ solar torque expression (2.14) 


reduces to 



= CcrCS-] - ^2) jsia(0+c+X+ 6^)) 


X [2 cos(e+5+X+ 6 ) cos 6 

0 © 

- Bin(e+c+X+ 6 ) sin 6 j (2.18) 

The simplified expression indicates that for small control 
surface rotations about the 'off position 6 , the radiation 
torque depends on the single control variable (6^ - 62). 

This fact can be used to advantage during mechanization of the 
system as the tw surfaces may be rotated in opposite directions 
by equal amounts employing a single drive arrangement. 

l*urthermore, the solar torque depends on the choice of 5 

© 

in addition to its time-varying character. Maximization of the 
torque producing ability of a given controller, characterized 
"by I Mg/C(&^ - ^2^ I » appears to be a suitable criterion for the 
selection of 5^. This however, leads to a 6^ dependent on both 
0 and X. On the other hand, it would be desirable to have a 
constant 6 about which the control surfaces rotate. Hence, 
maximization of the average torque generating ability over a 
torque cycle, 

I *2^ lav jsiu(0 + fig) {2 cos(0+6g) cos 6 ^ 

- sin(6+6g) sin 6^} { a:© (2.19) 

is considered, where 0 = (0 + 5+X) represents a dummy variable of 
integration. 
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It may be easily shown that the integral is a maxirmam for 

6 ^ = niT , minimum for 6 = ( 2 n+ 1 )ir /2 where n = 0 , 1 , 2 , .... 

6 0 

Taking 6 ^ = 0 , letting 6 ^ = 6 , 62 = “^ substituting 
Equation (2.18) in Equation (2.13), the equation of motion 
takes the fom. 

X " + 5K sin X cos X = 40*^5 | sin(6 +C4A ) 1 Gos(e+?+x) (2,20) 

The control law for the variable should be so chosen 
that the controller introduces positive damping and restoring 
effects in the pitch dynamics of the satellite. The control 
system should also exhibit asymptotic stability for small 
amplitude motion over a range of system parameters. To meet 
these objectives, we consider the control law 

d = - sgn { cos(9+?)} (yX*+vX) (2.21) 

The dynamics of the controlled system is described by the 
differential equation resulting from the substitution of 
Equation (2.21) in Equation (2.20). Eor pitch motion near the 
equilibrium attitude, linearization in 1 leads to 

X” + (ul sin 2 nj) X’ + (v {sin 2nj + 3 K)X = 0 (2.22) 

where y = 20 cy , v = 20 av and n = ( 6 + 5 ) is a new 
independent variable. 

It is apparent that the damping and restoring torques 
introduced by the solar controller are time-varying in 
character. However, they remain nonnegative throughout. One 



may, therefore, expect to find a range of combinations of the 
parameters p , v and K leading to asymptotically stable 
operation. 

2,5 System Stability and. Performance 

The time -varying coefficients in Equation (2.22) are 

periodic functions of p with a period of tr/2. Hence, Ploquet 
31 

theory may be applied to determine the stability of the 
system. 

letting x^ = X , X 2 = X ' » x = (x-j , X 2 ) , Equation (2.22) may 
be expressed in the form 


x’ = M(n ) X 


(2.23) 


where 


M(ti ) = 


1 


-(v| sin 2 n| + 3K) -p] sin 2 n| 


(2.24) 


The stability analysis involves the examination of the 
characteristic multipliers of the system (2.23). These are 
the eigenvalues of the square matrix Y(tt/2) where T(n) 
satisfies the matrix differential equation 


l' (n) = M(n) Y(n) (2.25) 

and the initial conditions Y( 0 ) = I, the unity matrix. The 
system is said to be as 3 rmptotically stable if, and only if, the 
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moduli of all characteristic multipliers are less than unity. 

It is unstable if any characteristic multiplier has a modulus 
greater than unity. 

In the absence of a closed -form solution, the system 
(2,23) was integrated numerically over the period 0 to 
n = •ir/2 to obtain the stability information. The results of the 
analysis were summarized in the form of stability charts in the 
y, V parameter space for different values of the satellite 
inertia parameter. The stability investigation extended over a 
large range of the controller gains. Por a clearer appreciation 
of the resulting stability-instability patterns, the stability 
charts for a restricted range of y , u are shown in figure 2.3 
for two extreme values of K. 

It is apparent from figure 2.3a that for a dumbbell 
configuration aligned nominally along the local vertical 
(K = 1), a wide range of parameters results in asymptotically 
stable motion. Instability occurs only for certain regions of 
small y values. These small instability regions were found to 
appear repeatedly in a pattern as larger values of v, beyond 
those shown in the figure, were considered. The most adverse 
situation of a dumbbell configuration nominally aligned along 
the local horizontal direction (£ = -1 ) is represented in 
figure 2.3b. The instability region is significantly enlarged 
as the gravity torque now acts as a destabilizing influence. 
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The appearance of instability near the y axis indicates that a 
minimtmi amo'unt of restoring torque from the controller (minimum v) 
is essential to have stable operation. Also, for small v values 
the instability region enlarges somewhat as y is increased, 
suggesting that an increase in the time-varying damping torque 
may at times lead to instability. Nevertheless, the results 
establish the capability of the solar controller in achieving 
pitch stabilization along the inherently unstable gravitational 
equilibrium attitude. The stability analysis for intermediate 
values of the inertia parameter showed similar stability- 
instability regions with the general trend of increasing stable 
regions as K increased from -1 to +1 in the physically possible 
range of -1 1 K 5 1 . 

Note that y and ^ would vary slowly as the time of the 
year changes. Their variation is relatively small for ' 
satellite orbits having a low inclination from the ecliptic 
plane. Hence, it would be desirable to select the gains y, v 
so that the system operating point lies well within the stable 
region of the stability chart. This would ensure asymptotically 
stable operation throughout the year. 

Another feature characterizing the controller's 
performance would be the speed of attitude correction it is 
capable of. The largest '’system time constant”, which provides 
an estimate of the time taken to reduce the amplitude or the 
amplitude envelope by a factor 1/e, may be taken as an 
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approximate measure of the controller’s speed of response. 

The largest time constant is given hy 

' 2 . 26 ) 

where is the largest of the moduli of the characteristic 

multipliers. The dependence of t on the controller gains and 
inertia parameter was found to be quite complicated. G-enerally, 

T decreased with increase in il, v and K. Time constants of 
the order of a few orbital degrees (« 15 - 20®) were found 
achievable within the parameter range 0 < y , v < 50. 

The full nonlinear equations of motion (2,13) were 
integrated in conjunction with the control law (2.21) for 
different values of the controller gains and initial conditions 
to study the system response. Typical responses obtained are 
presented in Pigures 2.4 and 2,5, These show the attitude 
responses and the control histories subsequent to an initial 
position error and an impulsive disturbance. The response is 
observed to be of a damped oscillatory character withy = 0,1 
and 1,0 (Figure 2,4). When the parameter is increased 
to 1,0, with V remaining the same, the response is found to be 
of an over damped nature (Figure 2,5). This characteristic of 
the present nonautonomous system appears to be quite similar 
to a second order system with constant coefficients. 
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The results exhibit time constants ranging from about 50® 

to 150 ® of the orbit. These, of course, are due to the very 

small values of the gains employed. Considerable reductions 

in the settling times were found possible with the use of 

increased M and v values. It may be noted here that the 

impulsive disturbance = 0.1 represents a much greater value 

32 

than is likely to occur due to micrometeorite impacts . 
Consequently, much smaller settling times would be expected in 
reality. The typical system response presented m^ be better 
appreciated with a numerical example. The solar parameter 

2 

value C - IO is achieved for a spacecraft having I = 500 Kg-m 

2 

in the geostationary orbit with a control surface area A = 1m 
and moment arms = 3m. 

2,4 Concluding Remarks 

The conclusions based on the present study may be 
summarized as follows ; 

(i) The analysis shows that it is possible to design a 
mechanically simple solar pressure control system for 
pitch control of earth-pointing satellites. The 
suggested control law for the control surface rotations 
is also very simple to realize. 

(ii) The controller is capable of stabilizing the spacecraft 
pitch attitude along gravitationally stable as well as 
unstable orientations. 
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(ill) Ihe control system assures asymptotically stable operation 
over a wide range of system parameters. 

(iv) The controller appears to be uselUl for applications wherein 
control times of the order of a few degrees of the orbit 
may be acceptable, such as, long-life scientific missions. 
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3. THREE- AXIS, ATTITUDE STABILIZATIOH OF 
eaeth-poihtihg spacecraft 

The previous analysis suggests that it should he possible 
to modify the controller geometry to generate solar control 
torques about three mutually perpendicular body-fixed axes. 

The two-surface controller* may then be utilized for three- 
dimensional attitude control of earth-pointing spacecraft. Th< 
development and analysis of a three-axis attitude control syste 
using only a single pair o£ rotatable control surfaces forms th 
objective of this chapter. The equations governing roll, yaw 
and pitch degrees of freedom are formulated for a triaxial 
spacecraft moving in a circular orbit of arbitrary inclination 
from the ecliptic plane. Optimal control theory is applied to 
synthesize a feedback control law for the differential control 
surface rotation which promises asymptotically stable syston 
operation. The optimal control law involves the implementation 
of periodic time-varying gains. Subopt imal control strategies 
promising considerable reduction of the associated computational 
requirements are explored. Numerical results are presented 
TA?hich establish the effectiveness of the control system in 
stabilizing the spacecraft attitude even -under severe gravity- 
gradieiqt conditions and external disturbances. 
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5.1 Pormulation of the Prohlem 

Pigure 5.1a shows the geometry of the orbital motion of 

an unsymmetrical satellite with its center of mass S moving in 

a circular orbit about the earth’s center 0. Ihe rotational 

motion of the spacecraft about its center of mass may be 

conveniently specified by a set of finite rotations indicated 

in figure 5.1b. The rotations perfoimed in the sequence, y 

about the z^-axis, B about the y^-axis and ^ about X 2 -axis, 

bring the orbital coordinates x , y , z into coincidence v;ith 

o o o 

the satellite principal axes x,y,z. The rotations are referred 
to as the roll (y), yaw (B) and pitch (X) angles, respectively. 

The proposed solar controller configuration is indicated 
in figure 5.2. It consists of two highly reflective lightvjeight 
surfaces P-] and P 2 with their center line P 2 located at an 
angle Xp from the y-axis in the yz-plane. The surfaces 
intercept the solar radiation pressure resulting in a torque ■ 
about the satellite center of mass. The solar radiation torque 
may be controlled by rotating the control surfaces about their 
body-fixed axes by angles 5^ and Bg S'S shown in the figure. 

The dynamical equations governing the attitude motion may 
be derived using the classical lagrangian approach. The' kinetic 
and the potential energies associated with the rotational motion 
of the spacecraft in the gravity-gradient field are found to be 

T = (1/2)(I + I + I 

% X ^ y y z. z-' 


(5.1) 




/ 

/ 



Fig. 3-2 Solar controller configuration 
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where 


• . • 

ti) = 6 cos Y cos 6 - Y sin 3 + X (3,2a) 

m 

mm ^ 

Wy = ($ - 0 sin y) cos X + (6 cos y sin 3 + 

Y cos 3 ) sin X (3.2'b) 

= -(3 - 6 sin Y ) sin X + (6 cos y sin 3 + 

t cos 3) cos X (5.2c) 


^ 0 o 

Y = -i^^4 R^) t (21^ - sin Y cos 3) 

+ { (sin^Y sin^ 3 - cos^ Y ) cos 2 X 

- (sin 2 y sin 3 ) sin 2X} 1 (3.3) 

The equations of motion are then obtained from 


(d/dt)(3L/3k) - (3I/3k) = Qj^ 

k = Y , 3 , X 


(3.4) 


where L = (U - Y) is the system Lagrangian and (k =y > 3 ) 

denote the generalized forces due to solar radiation pressure. 


The solar pressure torque about the satellite center of mass 
is given by 

M = e(cos X_ I + sin X_ k) X P. 

D y F ’ 

-e(Gos X T + sin X _ k) X (3.5) 

It . F 

where 

= -2P pAj cos j cos n^, i = 1,2 


(3.6) 



(3.7) 


n. = (cos + (sin ^ sin 

-(cos sin 5j^)k , i = 1,2 

Using the principle of virtual work, the generalized forces are 
obtained as ; 

Q = (2 P pAE) [ h^ cos 6 cos( X + X ) - ^2 

( 3 . 8 a) 

Q 0 = -(2ppAe) h 2 sin(^ + X^) (3,8h) 

= (2p pAe) h^ (3. Sc) 

Here, the quantities h^ and h 2 are defined to be ■ 

=| cos 5 -j I cos cos - jcos ^ 2 1 ?2 ^2 

(3.9a) 

h 2 = I cos Cj I cos ? ^ sin 6^ - j cos ^ 2 I ^ 2 ^ 2 

(3.9b) 

where 

cos cos + (e^ sin X^ - e^ cos Xp) sin 6 j_, 

(i = 1,2) (3.10) 

and the earth-sun unit vector components e^, e^, e^ are given 
in Appendix A, 

Substituting Equations (3.1), (3.3) and (3.8) in Equation (3.4), 
and nondimensionalizing the independent variable through the 
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relation 8 = Qt, the equations of attitude motion finally 
take the form : 

^ P sin^ X + cos^X)cos^ 6 ) y" 

+ {(I - I ) cos g sin X cos X} g” - (l sin 8) x" 

y ^ X 

= f^ + C {h^ cos g cos (X + X ) - h2 sin g} (3.11a) 

* Jr 

2 2 

{ (l„ - I„) cos gsin X cosx}y" ^ sin X}g'' 

y zt ' y z 

= f - Glx^ sin (X + X„) (3.11h) 

$ P 

- (I^ sin B)y”+ I^x" = ^X "^’^2 (3.11c) 

Here, the quantities fj^ (k = y , g , X ) are rather lengthy 
functions of the attitude angles and their derivatives. Their 
full forms are presented in Appendix A. The functions h-j , h2 
indicate the dependence of the radiation torque on the control 
surface rotations 6^ and 62- The parameter 0 characterizing 
the magnitude of the solar torque generated hy the controller 
is given hy 

C = 2pp Ae/ (3.12) 

The equations of motion consist of a sixth order, coupled, 
nonlinear, nonautonomous system of differential equations. 
Consequently, the selection of a suitable control policy for 
the control surface rotations 5^, Sg indeed formidable. 
Therefore, an approach based on the linearization of the 


i 



equations of motion about the system equilibrium configuration 
is presently considered. 


3.2 Control Synthesis 

An inspection of Equations (3.11) indicates that the 
system possesses the equilibrium configuration : 

Y = B = X = 0 

6^ = 62 = fig (arbitrary) (3,13b) 

Physically, this implies that v/hai the two control surfaces 
are parallel and no external disturbances act on the satellite, 
the body-fixed principal axes, x,y,z, would remain aligned with 
the orbital coordinates, x^, y^, z^, respectively. This 
corresponds to the nominal earth-pointing attitude of the 
spacecraft. The objective of the solar controller is to maintain 
the earth orientation of the satellite in presence of external 
disturbances. Eor small amplitude motion near the equilibrium 
position, the governing equations may be simplified through 
linearization. Substituting 6^ =5^+5^, dp " ^2 

A# 

linearizing in » Equations (3.11) take the form : 

y” - (Kg-I) B’ + 4K2Y = ( 0 /I^)g^(e,<|>, 6 ^)cos Xp(«i - 62^ 

(3.14a) 

6” - (K|+1) y’ - K-iB = ^g)sin Xp(i^-52^ 

(3.14b) 



■where 
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X" + 3 [ = 

(G/yg 2 (e, 4 », 5e)(V^2) (3.14c) 

(6 ,<|! , 6 ^) = ( 1 / 2 ) jb( I 3 a sin 26^ + 2 ' 

2 

(cos 6 ^ + cos 2^ ) sin (e +C )J (3.15a) 
0 0 

g 2 ( 8 ,<^, <S^) = ( 1 / 2 ) jh{ [ 2 a(sin^ 6 ^ - cos 2 6 ^) + 

3 ^/l-a sin 25 sin (9 +? )J ( 3 , 1513 ) 


with 


a = Sind) sin 1 


h = - a cos 6 + sin 5 sin(9 +c ) 

0 0 

C = X - tan” (tan <l> cos i) 

ir 


(3.16a) 

(3.16h) 

(3.16c) 


The linearized equations of motion reveal certain features 
of the dynamic structure of the system. Por small control 
surface excursions about the 'off' position, the radiation 
torque depends only on the difference ( 6 ^ - 62 ) between the two 
control surface rotations. Consequently, all the three degrees 
of freedom must he controlled by maneuvering the single variable 

mm 

( 6 ^ - 62 ). This may lead to somewhat complicated control 
synthesis. Prom a mechanizat ion point of view, however, it 
appears advantageous as pointed out earlier. 

Purthemore, the uncontrolled system exhibits a 
gyroscopically coupled roll-yaw (out -of-orbit -plane) motion 



and an uncoupled pitch (in-orbit -plane) motion. Hence, to 
achieve a good overall performance, the controller should 
provide significant torque contributions in both these modes. 
Hondimensional measures of the solar torques controlling the 
roll-yav and pitch motions per unit differential control 
surface rotation are given by 

“yB + Kg /C («, - = giCB.*. Bg) (3.17a) 

= ly^/G {? - Bg) = ggO,*, Sg) (3.17b) 

Por best perfomance, the 'off’ position 6^ of the control 



are as large as possible for all values of ® and 4> . Ihe 

optimum ^ would obviously be a function of both 6 and 4> , 

On the other hand, it v/ould be desirable from practical 

considerations to have a constant S about which the control 

e 

surfaces rotate. Acceptable constant values of may be 
obtained by requiring that the average values of both I MLygl 
and I I over an orbit be sufficiently high. Restricting 
attention to near ecliptic orbits (i < 30®), the effect of “J* 
on the magnitudes of g.j and g 2 becomes relatively small. 
Analysis shows that most choices of except very near 
= 0 or 90®, lead to sufficiently high values of the 
average out-of-plane and in-plane solar torque components. 

Best results were obtained with S s 30®, 105°. 
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Optimal control theory is next applied to derive a feedhaol 
control law for the differential control surface rotation 
This is rendered convenient hy first expressing Equations (3.14/ 
in the state variable form. letting x^ = y > ^2 = Y' , = 3 , 

x^ = g’ , x^ = X , Xg = X’ and n = they take the form, 

i’(e) = ibc( 0 ) + B( 6 ) u( 0 ) (3.18) 

where the matrices A and B are given by, 




J 



[ ?( 6 ) Q(0) xCe) + T? (e ) ¥(0) u(e)yd0 


(5.20) 


•where Q(9) and ¥(0) represent positive semidefinite and positive 

definite weighting matrices, respectively. The control law 
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minimizing the index (3.20) is well known to he 

u(0) = F(e) x(0) (3.21) 

where the gain matrix P(6) is given by 

P(0) = -¥"^(0) B^(0) K(0) (3.22) 

Here, K(0) represents the solution of the matrix-Riccati 
equation 

K’ (0 ) = -K(e ) A -A®K(0 ) - Q(0 ) 

+ K(e) B(0) ¥(0)“’' B®(0) K(0) (3.23) 

with boundary condition K(9^ "^ “ ) = 0. 

5.3 Computational Procedure 

The attitude control performance of the system was 
evaluated for a specific satellite equipped with the solar 
controller. The data used and the computational procedure 
employed are summarized below ; 

(i) A triaxial spacecraft in the geostationary orbit 

(i = 23.5°) with the following inertia properties is 
considered: = 400, = 275, = 620 kg - m^. 
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These inertia 'values lead to the shape parameters 

= 0.8, K 2 = 0.2 v/hich correspond to an unstable earth - 

35 

pointing configuration in absence of the controller 
The inertia parameter values were so selected purposefully 
in order to exspine the controller performance under 
adverse gravitational conditions. 

The control surface areas, moment arms, the body location 

of the control surface rotation axes and their 'off 

2 

position are taken tobe, A=1m,e = 3m,A^ = 0, 

5 = 105 ®, which result in the parameter 0 = 4860 kg-m . 

(ii) The choice of weighting matrices Q(®) and W(0), is somewhat 
arbitrary. Presently, they are both taken to be identity 
matrices. This corresponds to assigning equal weights to 
the attitude error, attitude rate and the differential 
control surface rotation. 

(iii) A value of the solar aspect angle <j> is selected next. 

This completes the parameter specification for the 
problem. 

(iv) The matrix-Riccati equation (3.23) is then solved backward 
in time with the boundary condition = 0. In 

practice, = 10 orbits was found to be sufficiently 

large for the solution K(0) to reach a steady state 
periodic character with a period , 



(v) The associated 2 it - periodic gain matrix, 

p(8) = iPy Py’ ^8 ^8’ ^X* ^ (5.24) 

is evaluated next from Equation. (3.22). Pourier 
coefficients are fitted to the system feedback: gains to 
facilitate their subsequent generation. 

(vi) Pinally, the full nonlinear Equations (3.11) governing the 
system dynamics are integrated to obtain the attitude 
response of the system. The simulation requires the 
control surface rotations 8^, ^2 measured from the yz-plane. 
Assuming the surfaces to be driven by equal amounts in 
opposite directions from the 'off position 5 one obtains, 

6^=6^ + u/2 (3.25a) 

(3.25b) 

The possibility of loss of control during the eclipse of 
the satellite by the earth was also considered. Por this 
purpose, the shadow criterion (Appendix C) was introduced 
during the evaluation of the attitude response. 

(vii) Por different values of 4> , the steps (iv) to (vi) are 

repeated to evaluate the system performance at different 
times of the year. It is noted that replacing <j> by 

+ ^) and u by - u leaves the governing Equations (3.11) 
unchanged. Hence, it is only necessary to vaiy over a 
half-year period. 
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3.4 Results and Disoussion 

figures 3.3 and 3.4 present the system gains for different 
values of the solar aspect angle , The gains are plotted as 
functions of h = (6 + c) for convenience. It is apparent that 
the nature of the pitch gains is quite different from 

that of roll and yaw gains F f , and Po ,P«». In 

Y 7 Y p ^ p 

particular, Pg always maintain opposite signs while 

the corresponding rate gains P^i and Pgt always have the same 
sign. This hehaviour appears to he due to the gyroscopic 
coupling between the roll and yaw degrees of freedom. 

Although the gains has^e been computed for the satellite 
in the equatorial orbit, the gain histories shown for <{> = 0 
have a much wider implication. Mote that the governing system 
Equation (3.'i4), with n as the independent variable, remain 
the same in the situations i = arbitrary, = 0 and i = 0, 

<f> = arbitrary. Consequently, the results shown for <P = 0 

also represent the gain histories required for a satellite in 
the ecliptic plane at all times of the year. In such a case, 
the gain histories show a skew -symmetric variation from their 
zero mean value which corresponds to the skew- symmetric 
nature of the functions ) , §2^®^ with the parameter 
a = sin 4> sin i equalling zero. 

As pointed out earlier Pourier series representat ion for 
the periodic gains was employed for evaluating the system 
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variations of 4> in the latter range of <1> (figures 3.3 and 3.4). 
The analysis thus suggests that a judiciously selected gain 
scheduling scheme may he introduced to achieve suhoptimal 
control resulting in reduced gain computations. 

In the present analysis, the yaw control is achieved 
entirely through the coupling between the roll and yaw motions 
as the controller location X was taken to he zero. ¥ith a 

ir 

suitable selection of X , a better distribution of the control 

ir 

action betv/een the roll and yaw motions v;ould be possible with 
the prospect of an improved overall performance. In general, 
the control surface rotation axis may be skewed at any angle 
relative to the spacecraft axes x,y,z. Its best orientation 
would depend on the orbital inclination from the ecliptic. 

3.5 Conclud in g Remarks 

The significant conclusions based on the analysis may be 
summarized as follows ; 

(i) A solar pressure attitude control system for the simultaneous 
roll, yaw and pitch control of earth-pointing spacecraft has 
been developed. The controller design minimizes the hardware 
requirements as it involves only two control surfaces which 
may be rotated differentially by a common drive mechanism. 



54 


(ii) The hardware reduction is associated with the implementation 

of time-vaiying periodic gains which need to he continuously 

updated as the sun -line rotates during the year. However, 
effective suhoptimal control may he achieved hy adopting 

a judicious gain scheduling scheme which involves updating 

the gains only a few times per year. 

(iii) The controller, even with moderately sized solar surfaces, 
is found to he capable of maintaining the spacecraft 
attitude again^ severe external disturbances as well as 
destabilizing gravitational toroues. 

(iv) The control system promises asymptotically stable operation 
throughout the year. 

(v) The influence of the earth's shadow on the performance is 
found to be negligible for spacecraft in geostationary 
orbits. 

(vi) The perfomance of the system is likely to improve through 
optimal orientation of the control surface location 
relative to the spacecraft body. 
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4. PIE OH STABILIZATIOF OP SPAOECRAPI IF M 
IIERTIALLY-PIXia) AEMUPE ‘ 


The analyses so far have dealt with the prohlem of 
achieving an earth-pointing spacecraft attitude hy means of 
solar radiation pressure. However, its utilization for attaining 
a space-oriented satellite attitude has remained unexplored. On 
the other hand, there are long-life spacecraft missions, such as 
those involving astronomical ohservations, which demand stahili- 
zation along inert! ally -fixed orientations . It therefore 
appears interesting to extend the solar pressure control tech- 
nique for inertial attitude stabilization. If successful, it 
would represent an inexpensive method of attitude control for 
an important class of space applications. 

This chapter examines the potential of the two-surface 
solar controller for stabilizing the pitch attitude of an 
unsymmetrical satellite along an inertially-fixed orientation. 

A control strategy involving a nominal control to counter the 
gravitational torque and a feedback control to impart stability 
to the system is developed. An essentially analytical approach 
leads to control laws that are easy to implement. The stability 
properties of the system are established for all times of the 
year. The influence of the important system parameters and the 
earth's shadow on the system perfomnance is analyzed. 
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4.1 Problem Pormulation and Control Synthesis 

Pigure 4.1 shows the geometry of the orbital and attitude 
motion of an unsymmetrical satellite with its center of mass S 
moving in a circular orbit about the earth's center 0. The 
inertial coordinate system X,Y,Z is selected such that the I- axis 
lies at an angle a from the line of nodes, YZ represents 
the plane of the orbit and the X-axis -is along the orbit normal. 
The body-fixed principal axis y points in Y direction when the 
spacecraft has the desired in ertially- fixed pitch orientation. 

The orbital angle ® is measured from the Y-axis. The angle a 
represents the inertial pitch attitude of the spacecraft measured 
from the line of nodes. The controller configuration is assumed 
to be the same as in Chapter 2 for maintaining an earth-pointing 
attitude. 

The equation of pitch motion may be readily obtained by 

replacing 6 by ( 8 + a.^) and A by ( a - a - 6 ) in Equation 

© © 

(2,13) as 

ct" + 3K sin (a - a - 6) cos (a •- a - 8) = 

© © 

c sin(^ + a- sin(C + a _ + 5^) cos 

- ( sinCc + a - Cg + ^2^ I sin(C + a - a ^ + 6^) cos 62} (4.1 ) 

where the solar parameter C, characterizing the controller size, 
is defined as 

C = 2P pAe/fl^I^ 


(4.2a) 
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and the quantities a and ? are given hy 

a = 1 - sin% sin^i (4.2h) 

?C<J>) = - arotan (tan<|) cos i) (4.2c) 

In absence of the controller, a = a + 6 satisfies Equation 
(4.’^). Physically, this corresponds to gravitational stabiliza- 
tion in an earth-pointing attitude. Inertial stabilization in a 
fixed attitude would require a =ag to bean identical solution 
of the pitch equation. Setting a = a in Equation (4.1), the 
nominal controls required to coimter the gravity torque are 
governed by 

|sin(C + ^-^q) I sin(C + 6^^) cos 6^^ - jsin (5 + ^20^^ 

sin(C + 5 20^ ^20 ~ -(5K/2Cff)sin 20 (4.3) 

The two unknowns 52 o( 8), however, cannot be determined 

from this single relation and a degree of arbitrariness remains 
as to their choice. A convenient constraint would be to permit 
them equal and opposite rotations about some mean position 6^, 
i. e. , let 

= 5g + «„(9) (4-4a) 

«2o(e) = «e - (4.4b) 

Substitutiag in Equation (4.3), the nominal control S^^CS) Is 


then obtained from 
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I sin( C + <5 + 5 )1 sin(C + 5 +. 5 ) cos(5 + 5 ) 

fc* O © 0 © o 

- |sin( c + 5 - 5 )| sin( C + 6 - 6 ) cos( 6 ^ - 5^) 

O * 6 0 6 0 

= -(3K/2Ga)sin 26 ( 4 . 5 ) 

The transcendental nature of Equation (4.5) indicates that, 

in general, continuous onboard numerical computation must be 

performed to solve, for 6 ^( 0 ). On the other handj an analytical 

solution -would be much easier to implement. It is apparent that 

5 ^ = 0 represents the solution -when the right hand side of 

Equation (4.5) is zero, as wo-uld be the case for a satellite vdth 

symmetry about the pitch axis. When conditions exist such that 

1 throughout the year, the solution 5 ^( 6 ) ■'would be 

expected to remain small. Noting that 1 k 1 1 and ^ > cos i, 

the requirement can be met for most satellites -with a fairly 

moderate controller size, except in orbits -with "too large an 

inclination from the ecliptic. As an example, for a satellite 

2 

in the geostationary orbit having I = 250 Kg-m , K = 0,5 and a 

o 

controller with A = O .25 m , a = 5m, = 0.114'^< 1. 

Eor small 5^(6), therefore, the left hand side of Equation (4.5) 
may be expanded in ^o\-}ev8 of 6 ^ and only the linear terms 
retained. This results in the approximate solution 

6 ^( 6 ) = -(3E/4CaD) sin 26 ( 4 . 6 ) 

where 

E = I sin(? + H 2 cos 6 ^ cos(? + 5^) 

- sin 5^ sin(C + ^ 

6 , , 6 


(4.7) 
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The analytical approach ■would he justified pro'vided the 

Quantity D does not hecome too small at any time of the year. 

The obvious solution lies in suitably updating the mean 

position 5 as C varies slowly, i.e. , let 6 =6 (C). Best 

results would be obtained by upgrading 6 such that it maximizes 

} dI for a given The variation of the optimal § with c is 

0 

indicated in Figure 4.2a. The relationship is not only periodic 
with a half-year period but is also practically linear, making 
the selection of the optimum 6^ quite simple. The figure also 
shows the variation of the associated jD | which remains within 
the limits '1.16. This establishes the validity of 

the analytical solution for all times. Its accuracy was assessed 
by comparison with the exact solution obtained by numerically 
solving Equation (4.5) for several parameter combinations. 
Excellent agreement v;as found as reflected by the typical plots 
shown in Figure 4.2b. 

It should be pointed out here that various other choices 

of the relation 6 = d (C ) would also result in sufficiently 

0 0 

large values of |d|. For example, 5^ may be assigned a 
constant value over a range of C(i.e., ) and then updated to 

a new constant value over the next range of ?. This would 
amount to replacing the optimum cuive in Figure 4.2a by a 
stepped function. 

Although the nominal control opposes the disturbing 
gravitational torque, it does not stabilize the satellite 
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against initial errors or disturbances such as might arise from 
micrometeorite impacts* A suitable feedback control system is 
hence required, letting 6 denote the additional feedback 
controlled rotations to be provided, the total control surface 
rotations become 

+ <5^ + 6 (4.8a) 

‘ o 

dp = «e " '’^0 " ^ (4.Sb) 

Por the general case of an arbitrarily large 6 q( 9)» the 
equation of motion in the neighbourhood of the desired 
inert ially-fixed orientation a may be obtained by letting 
a = + a , substituting Equations (4.8) in Equation (4.1), 

and linearizing in a and 6 . Subsequent to considerable 
algebraic manipulations, with due regard to Equation (4.5), 
this results in 

a” + {5K cos 26 - Caf^(0)}S - {Ca f 2 (e )}6 i 0 (4.9) 

where 

f , (0 ) = 2 {| siniz + d + 6 )| cos(C + d + 5 ) cos( d + 6 ) 

-I sin(? + dg - 6 ^) j cos(c + " *^ 0 ^ oos(S^ - d^)} 

( 4 . 10 a) 

f 2(6 ) = f^(e) -•flsin(? + 6 g +'5o)l sin(c + d^ + 5^)sin(fig+5^) 

-1 sin(c + dg - 5^)1 sin(? + d^ - 6 ^) sin(5g - 5^)) 

(4.10b) 
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station (4.9) represents a linear controlled system with 
oomplxoated time-varying coellloients dependent „n the nominal 
control d„(a). low order ol the system, however, suggests 

the possibility of analytically synthesising a suitable feed- 
back control law. 


■ The problem simplifies considerably if attention is 

restricted to the case of small E/Oavalues leading to small 

d„(9). one may then expand the coefficients of J and d m 

Equation (4.9) and retain only the first order terns In d to 
obt ain ° 


t {3K oos 26 . 20a Ed^); . ( 20 a 1)4 , q 

where 

E={2 cos d^ oos 2(c. d^) . ^ 

sgn { sin( c+ 5 ^)} ( 4 . 12 ) 

Substituting the analytical solution (4.6) for <S ( 0 )^ the 
goV6rning ©QUci'fcion ■fch© fopm 

5" +{ 3E(1+E^/4d 2)^ cos(2e - v)}a - (2CaD)6 . o ( 4 .I 3 ) 
where 

V = arctan (E/2D) 

(4.14) 

In absence of the feedback control 5, the attitude motion 
near the nominal orientation is governed by an undanped second 
order differential equation with a periodic coefficient. The 
control law for S should be so chosen that it introduces 
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corrective restoring and damping effects into the system. A 
logical choice ■which is simple to realize woiild he a piecev/ise 
linear control law of the form 

6 =-(ua* +va) sgn D (4.15) 

where \i , v represent constant system gains. 

4.2 Stability Analysis 

The usefulness of the proposed control law would he 
established pro'vided it leads to asymptotically stable system 
operation for a wide choice of the controller gains. 
Substituting Equation (4.15) in Equation (4.13), and shifting 
the independent variable from 6 to ( e - V2) for convenience, 
the dynamics of the controlled pitch motion is governed by 


a"+ya’ 4-(v + q cos 20)5 = 0 (4.16) 

where the parameters u , v and q are defined as : 

U = 2Gor !d1u (4.17a) 

V = 2C ojll V (4.l7h) 

q = 3K;(1+E^/4I)^)^ (4.17c) 


Equation (4.16) represents a linear differential equation 

51 

with a time-varying periodic coefficient. Eloquet theory 
may therefore be applied to investigate the stability 

'Mr 

properties as a function of the parameters U and q. The 
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method involves the examination of the' characteristic 
multipliers of the system. In absence of a closed foim solution, 
their determination requires numerical integration of Equation 
(4.16) t'wice over the period tr for every combination of u , v 
and q. Considerable saving of computational work occurs if one 
eliminates the first derivative term through the transfomation 
~ - y0/2 

a = X e ' (4.18) 

which converts Equation (4.16) to 

x" + (h + q cos 26 )x = 0 ( 4 . 19 ) 

where 

h = V - y^/4 (4.20) 

Equation (4.19) is the classical Mathieu equation whose 
stability properties are well known. In the h,q plane, the 
stable and unstable regions occur alternately with transition 
curves associated with periodic solutions of period it and 2ir 
separating them. Eor a chosen value of q, a sequence of values 
of h separating the stable and unstable regions is obtained. 
These in turn map into parabolas in the y, v plane (Eigures 4.3 
and 4.4) in accordance with Equation (4.20). i-s the solution x 
can at best exhibit ordinary stability, the asymptotic stability 
of a requires that y, and hence the gain y , must be positive. 
The entire region of stability for x thus corresponds to 
asymptotic stability of a. 
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To determine the regions of asymptotic stability of a 
which o\7erlap the instability regions of x, consider the 
solution X of Equation (4.19) which has the form 

Yi© , Yp0 

X = C^.e p^ ( 0 ) + 02© P2(0) ( 4 . 21 ) 

where p-^(6), periodic functions of period ■jr, 0^ and 

Gp represait arbitrary constants and Y^ , Yg are the 
characteristic exponents of the system. The real parts of 
Y-[, Y 2 obtained from the characteristic multipliers 

Pl f p2 

r^ = In (p^ j /it , i = 1 , 2 (4.22) 

Substituting Equation (4.21) in Equation (4.1S), it is apparent 

that in the instability region for x, where at least one of the 

r. > 0, a 0 as 9 « provided u/2 exceeds r^^^, the 

1 max 

greater of the r^. 

With the selected value of q, a number of values of h 
lying in the unstable regions of x in the h,q plane were 
identified. In conjunction with numerical integration of 
Equation (4.19), Ploquet theory was used to determine the 
characteristic multipliers and the real parts of the 

^ -MW 

characteristic exponents. Proceeding to the y , v plane, this 
results in parabolic contours of constant filling the 

instability regions of x. The asymptotic stability region 
for a , overlapping the instability region of x, is then 
identified by requiring that 11/2 > 
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Ths slj slDilitiy invGs’fcigstiou was cairrisd. ouij foi^ s©ve 3 ?aX 
values of q and the resulting information is presented in 
Figures 4.3 and 4.4 in the form of stability charts in the 
il, V plane with q as a parameter. The stability charts may be 
thought of as cross-sections of a three-dimensional y , v , q 
space depicting the stability behaviour. Fote also that, in 
view of the symmetiy properties of the Mathieu equation, the 
same charts result for both positive. and negative values of q. 
As anticipated, the instability regions increase with an 
increase in q. 

Proper interpretation of the stability results requires a 
consideration of Equation (4.1?) which indicate that the 
parameters P , v and q acquire different values at different 
times of the year. Also, the nature of their slow variation 
depends on the upgrading policy 6=6 (?) through the 
quantities E and E. For a given satellite and controller 
(constant K, 0,u,v), the system representative point therefore 
moves during the year in the P, v plane and the stability- 
instability regions also change as q varies. A safe design 
would be achieved when the controller parameters C,y,v are 
such that the representative point remains within the stable 
region which is common for all the q values encountered during 
the ,|ear. Fortunately, the stability-instability pattern 
vaT'-'i4s with q in such a manner that the desired common region 
Ciir|H)5 readily identified by considering the stability chart 
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for the worst time of the year corresponding to the maximum 
value of q. 

For example, in cononnction with the strategy 6=6 ( 5 ) 
developed presently (Figure 4 . 2 a), a loose hound on the traverse 
of the repr esaitative point is easily obtained as the rectangle 
defined by 

20 li cos^i< ii < 2-32 Op (4.23a) 

2G Vcos\< V < 2.32 0 V (4.23b) 

Also, the greatest q value is found to be 

= 3-034 K (4.24) 

For a satellite representing the most adverse situation of the 
greatest gravity -gradient torque, }k{ = 1 which gives 
Oonsequently , any selection of 0 , p and v such that the 
rectangle described by Equation (4.23) lies in the stable 
region of Figure 4.4b would lead to stable operation throughout 
the year, A typical bound for a satellite in the equatorial 
orbit (i = 25 . 50 ) -with 0 = 5 , 0 . 25 , v = 0.6 is indicated 

in Figure 4.4b. It is apparent that a wide choice of these 
parameters is possible even under such conditions. With a 
reduction in K, the value reduces and the allowable 

range of the design parameters increases further. 

4.3 System Performance 

-cki estimate of the pitch damping rates the controller is 
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capable of may be obtained from the system 'time constant’. 
In the region of stable x, the envelope of a would decay 
exponentially with the time constant 


T = 2/5 = 1/yc|D|. (4.25) 

It is interesting to note that the maximization of { D 1 , 
which earlier led to a smaller nominal control 6^(9), results 
in quicker damping of the transients. The damping rates also 
improve with increased controller size and velocity gain. 


When the system operates at a point lying in the region 
of unstable x, the time constant increases somewhat to 

T = = VCuOlBl- (4.26) 

In any of these regions, the greatest value of r equals 

nisx 

half the height of the instability region for a, which is 
fairly small and reduces further with an increase in v. 
Therefore, in general, the increase in the time constant would 
be quite insignificant for larger V and/or v . 


The attitude response was evaluated by numerically 
integrating the nonlinear pitch equation (4.1) in conjunction 
with the control laws (4.6 and 4.15) and the shadow criterion 
(Appendix C) for various combinations of the system parameters. 
Initial conditions representing initial misalignments as well as 
the effect of impulsive disturbances were simulated. In orbits 
lying entirely outside the earth's shadov/, perfect inertial 
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stabilization was achieved subsequent to the damping of the 
initial disturbance in a fraction of an orbit. When the 
satellite encounters the earth's shadow during its orbit, a 
small attitude deviation following each passage through the 
shadow region was observed. Ihe magnitude and duration of the 
error were generally very small, and were found to be strongly 
dependent on the magnitude of the gravity torque and the length 
of the shadow period. 

figure' 4.5 presents the system response under the most 
severe earth's shadow condition which corresponds to the 
coincidence of the sun-line with the line of nodes ( <!> = 0, it ). 
The associated control histories are also indicated. Ihe 

parameters represent a satellite in the geostationary orbit 

2 2 
with I = 250 Kg-m and the controller dimensions A = 0,25 m , 

e = 3 m. It is apparent that the controller effectively removes 

the initial pointing error of 20® in less than half an orbit. 

The earth' s shadow effect leads to a periodic pointing error 

arising once per orbit. The maximum attitude deviation is found 

to be about 1 ® for the case of K = 0,5 (figure 4 . 5 a) . When the 

worst possible combination of the greatest gravity torque 

(K = 1 ) and the longest shadow period is considered, the error 

amplitude increases to about 2® (figure 4. 5b). The controller, 

even under such adverse conditions, is able to reorient the 

spacecraft quite quickly on emergence from the shadow. The 
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Fig. 4*5 Typical systtm response (a) K =0*5 
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performance may be further improved through an increase in the 
controller size and an optimum selection of the system gains. 

It may also be pointed out here that most orbits at the 
geostationary altitude would either be entirely free from the 
earth's shadow or encounter relatively small shadow periods. 
Consequently, the shadow induced error would be significant 
only during a small fraction of the satellite’s life-time. 

4.4 Concluding Remarks 

The . conclusions based on the present study may be 
summarized as follows ; 

(i) The feasibility of attaining an inert ially -fixed pitch 
attitude of a high altitude satellite by means of solar 
radiation pressure has been d ©non st rated, 

(ii) The analytical solution developed for the nominal control 
is very accurate and renders the control law simple to 
realize. The approximate solution remains valid for a 
variety of the updating strategies for the mean position 5^. 

(iii) Asymptotically stable behaviour of the system throughout 
the year over a wide range of the system gains is 
particularly attractive from the desiga point of view. 

(iv) The system performance, even under the most severe 
gravitational torque and earth shadow conditions, appears 
to be acceptable for a variety of missions. It may be 
further improved through an optimal selection of the 
controller gains. 
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5. STABILIZiTIOU OP SPACEGRAPP SYMMETRY AXIS 
ALONG AN INERPIALLY-PIXED APT I PULE 

, The previous investigation clearly establishes the 
feasibility of inert ially -fixed pitch stabilization by means 
of the two -surf ace solar controller. Phis chapter extends the 
approach to the problem of imparting any inertially-fixed 
orientation to the symmetry axis of an axisymmetric spacecraft, 
Por this more complex problem, a numerical approach is followed 
throughout. Phe two-dimensional character of the system 
necessitates the use of independent nominal rotations of the 
two control surfaces. Their differential rotation, however, 
may still be employed as the feedback control. An optimal 
linear feedback control strategy is synthesized through the 
minimization of a quadratic performance index which assures 
asymptotic stability for the system. Phe possibility of 
simplifying the software implementation through the use of 
subopt imal control policies is also explored. Phe effectiveness 
of the control system in stabilizing the spacecraft along 
arbitrary inertial orientations is evaluated through simulations 
of the governing nonlinear equations of motion. 

5»1 Pormulation of the Problem 

Pigure 5.1a. shows the geometry of the orbital motion of an 
axisymmetric satellite with its center of mass S moving in a 
circular orbit about the earth’s center 0. x^, y^^, z^^ represent 
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Fig- 5-1 (a) Geometry of orbital motion 
(b) Geometry of attitude motion 
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a system of coordinates parallel to the inertial frame Z,Y,Z 
hut with its origin at S, Two secjiential finite rotations, ot ^ 
about the z^-axis and otg about the y^ -axis are employed to 
bring the x^, y^, z^ frame in coincidence with the principal 
axis system x,y,z. Relative to the x,y,z coordinates, the 
'spacecraft spins about its symmetiy axis (x-axis) with an 
angular velocity The angles and Wg completely 

specify the attitude of the satellite symmetry axis in inertial 
space. Pigure 5.2 shows the two-surface solar controller 
configuration with the control surfaces located along the 
symmetiy axis of the spacecraft. The control surfaces are 
permitted rotations and measured from the xz-plane, 

about the x-axis. 

Presently, the Eulerian formulation is used to derive the 
equations of attitude motion in the gravity -gradient field. 

The Euler dynamical equations for the axisymmetric spacecraft 
may be writt en as ; 


to) X 


11 

5) y ” 0 

- I) (1) w 
^ Z X 

II 

“z + (1 

3 

3 

M 

! 

= (M 

gz 


SX^' X 

^sz>/"z 


(5.1a) 

(5.1b) 

(5.1c) 


where M , , M (k = x,y,z) represent the components of the 

gKl SxC 

gravitational and radiation torques along x,y,z axes, 
respectively. The angular velocity components are easily 
expressed as ; 
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(I) 


X 


- ct^ sin 02 + 


“y = 


a> 


o 


o^ cos 02 


(5.2a) 

(5.2b) 

(5.2e) 


The derivation of the gravitational torQue coinponents is 
quite lengthy but straightforward. With the usual first 
approximation including terms upto the order of (satellite 
dimension/R)^, they are given by^^j 


where 


M 

gx 

= 0 

(5.3a) 

gy 

= - 3 s?^(1 - I) I r r 

^ y X z 

(5.3b) 

gz 

= 3 - I) I r r 

' J -x J 

(5.3c) 

r 

X 

= (sin o^ cos 02 ) cosQ - (sin a 2 )sin 6 

(5.4a) 

r 

y 

= (cos a^) cose 

(5.4b) 


= (sin o^ sin 02 ) cos® + (cos 012 )®^^ ® 

(5'. 4c) 


Eollowing the procedure adopted in earlier chapters, the 
solar control torque components are easily found to be i 

(5.5a) 

(5.5b) 

(5.5c) 


M = 0 

sx 


M 


sy 


2 pp A e h^( « 2 ) 


M 


sz 


“ ^ P P ^ ^ ^ ^ ^ ^ 2 ^ 
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where 


h^(5^,52^ = I 

I cos 


cos 

^1 

cos 

9 


- 

I cos 

E 2 I 

cos 

52 

cos 

®2 

(5.6a) 


I cos 

«1 

I cos 

5l 

sin 

<^1 



jcos 

E2l 

cos 

52 

sin 

62 

(5.6b) 


Here, the angles of incidence are given hy 

cos = - e^ sin 6^ + e^ cos 6^, (i = 1,2) (5.7) 

and the satellite - sun unit vector components take the form : 


®x “ ^ (sin i cos cos a2 + cos i sin a2) 

+ cos <|) (sin a^j cos 02) 

Sy = sin 4) (sin i sin ot^) + cos 41 (cos a-j) 

= -• sin 4 ' (sin i cos sin ot2 ~ cos i cos 0 ^ 2 '^ 
+ cos4» (sin sin 02) 


( 5 .Sa) 

( 5 .Sh) 


(5.8c) 


Substituting Equations ( 5 . 2 ), ( 5 . 3 ) and ( 5 . 5 ) in the 
Euler dynamical relations ( 5 . 1 ), the desired equations governing 


a^, ap snd degrees of freedom are obtained as : 


■^(- Oj sin a2 + a^) = 0 


( 5 . 9 a) 


02 - (1 - I) Oj cos a2(- £*4 sin a2 + 05) = 

= _3 j2^(1 - I)r^ r^ + (2pp Ae/I^)h^ (fi-j , ^ 2 ) (5-9b) 



81 


*a^ cos 02 + *ai “2 + (1-1) (_ sin Og + ^) 

= 5 fl^(1-I)r^ Ty + (2p p Ae/Iy)h2(6^, 6^) (5.9c) 

Since ihe tortjue coinponenijs in ihe o,^ degree of freedom, 
are zero, it represents a cyclic coordinate. The associated 
first integral of motion given hy 

-aJ| sin = k (constant) (5.10) 

may he used to eliminate ct^ from the remaining equations of 
motion. The parameter k represents a nondimensional measure of 
the nominal spin rate of the spacecraft about its symmetry axis. 
For a nonspinning satellite, however, it would be zero. The 
equations describing the attitude motion of the symmetry axis 
finally take the nondimensional form ; 

a !j cos <*2 - a]| sin 02 + ("*-1)^ ap ■* 5(1-I)r^ r^ 

= G hgCS^, 62) (5. 11a) 

- (1-I)k alj cos 012 + ~ ^ ^2^ (5.11b) 

where the solar parameter G is defined as 

G = 2 p p Ae /I (5.12) 

The governing equations represent a fourth order, 
coupled, .nonlinear, nonautonomous, system of differential 
equations. The manner in which the control surface rotations 
and ^2 appear in the equations of motion is also quite 
complicated which makes the control synthesis problem quite 


difficult. 
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5.2 Control Synthesis 

It Is apparent from the equations of motion (5.11) that 
in absence of the controller and any external aistnrbanoes, 
ftj = aj = 0 represents an equilibrium configuration for the 
satellite. This, however, represents a specific case of an 
inertially-fixed symmetiy axis attitude with the symmetiy axis 
nomal to the orbit plane. An arbitraiy inertially-fixed 
orientation, a., = a,^, in general, would correspond 

to a nonequilibrium orientation in the gravity-gradient field. 

If the solar surfaces are so operated that the radiation torque 
exactly balances the gravitational torques, it should be 
possible to stabilize the symmetry axis along any orientation 
in space* The required nominal controls 6 .jq( 9 ) and 620 (e) would 
have to be such that = “pe ^®P^® 5 ents an identical 

solution of Equation (5.11). This leads to : 


hl(^lO» *^20^ = I 5(1-I)/G] r^g(e) r^^(e) (5.13a) 

^20^ = -[3(1-I)/C] r^^(9) r^g(9) (5.13b) 

where 

^xe^®^ “le ® (5.14a) 

r^g(e) = (cos a,,g)cos e (5.14b) 

^ 2 q(®) = (sin a.j^ sin a2g)cos9 + (cos a 2 g)sin 9 (5.14c) 


Ihese two nonlinear algebraic equations may be solved for the 
nominal controls 6-jo(6) and '^iew of the periodic 
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nature of the right hand +v 

s nana side, the nominal controls would be 

periodic functions of period tr . 

In addition, the solar controller would have to stabilize 

the spacecraft against disturbances such as due to Initial 

misalignments and miorometeorlte impacts.' Incorporation of 

suitable feedback controls jn +■: + 

xoxs, in addition to nominal controls 

would therefore be required. p„r small amplitude motion near 

the nominal orientation, the governing nonlinear equations may 

be simplified through linearization. Substituting 


“1 ®l e * “2 


-“2e -"“c. « 


2' “1 ' ho + h> *2 - ho + «2> 


^ It- C.KJ 

and linearizing In a, , Equations ( 5 . 11 ) 

become : 


ot'l + + a^a^ + 6., + a^ Sg = 0 

a^ + b.^S- ,.b2^^ -.b^a^ . b^ +b5 «2 = 0 


( 5 . 15 a) 

( 5 . 15 b) 


Here, a^, b^ (i = 1 , 5 ) represent very lengthy functions of 

^ ) <t> t **1e* *^2e’ ^10^^^* ^20^®^ system parameters. 

Their full forms are presented in Appendix B. 


Although independent control laws governing and 
be synthesized, it would not be necessaiy for attitude 
stabilization in view of the coupling between the a.j and S2 
iJiotions. On the other hand, speei^ing a suitable constraint 
between 6^ and 62 would reduce the problem to one with a single 
Control variable leading to considerable software simplifl.cat ions. 
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A careful examination of the functions a^ and a^ provides 
significant guidance as to the choice of the constraint. 
Equation ( 5 . 13 ) reveals that there exist values of e near which 
and 620 nearly equal irrespective of the values of I 
and G. Whenever the controller size is large (large G) or the 
spacecraft is near spherical (l = 1), it indicates the 


nominal controls and 62^ to he quite close to each othar 
for all 0. Equations (B- 4 , B- 5 , B -9 and B-IO) show that a^ 
and a^ (and similarly h^ and h^) would then tend to become 
negative of each other. Gonsequently , to generate larger 
feedback control torques, 6^ and 62 should have opposite signs. 
Bor convenience, they are presently constrained to be of the 
same magnitude also, i.e., 6^ = - 

Letting x.| = , ^2 = a-|> x^ = a2, x^ = , 

u = 5 ^ = - 6 ^, Equations ( 5 . 15 ) may be expressed in the state 
variable form as 


X’ ( 0 ) = A( 0 ) x( 0 ) + B(0) u(0) ( 5 . 16 ) 

where A and B matrices are given by 


” 0 

1 

0 

0 


0 


0 

— a^ 

-a^ 


^5 " ^4 

0 

0 

0 

1 

and B(0) = 

0 

-b2 

-b-( 

1 

-L3 

0 


?5 - "^ 4 . 


A(0) = 


(5.17) 
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Note that the elements of A and B matrices are periodic 
functions of the orbital time 0 -with the period tt . 

In order to obtain a feedback control law which assures 
asymptotically stable operation^-" of the system, we consider 
the minimization of the Quadratic performance index 

ef -»• «> 

J = (1/2) t3?(0) Q(0) x{Q) + i?(e) w(e) n(e)] d6 (5.1s) 

0 

where Q(0) and ¥(0) represent positive semidefinite and positive 
definite weighting matrices, respectively. The control law 
minimizing this performance index is well known to be^^ 

u(0) = B(0) i(0) (5.19) 

with the gain matrix I'(0) given by 

^(0) = -¥"^(0) B^(0) K(0) (5.20) 

Here, K(0) represents the solution of the matrix - Riccati 
equation 

K' (0) = -K(0) A(0) - A^(0) K(9) ~ Q(0) 

+ K(0) B(0) ¥*■^( 0) B®(0) K(0) (5.21) 

with the boundary condition = 0. 

5.5 Computational Procedure 

The performance of the proposed control system is studied 
through a response analysis of the system. A specific 
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satellite and controller configuration V7ith the following 
specifications is selected for the computations ; 

= 167 Kg-m^ 1^ = 1^ = 500 Kg-m^ 

A=1m^ e=3m i = 25.5'‘ 

This leads to the inertia parameter la 1/3 and the solar 
parameter C - 10 for a geostationary orhit . The computational 
procedure may he summarized as follows : 

(i) The appropriate value for the solar aspect angle <f> , 
determining the time of the year, is selected. 

(ii) Equations (5.15) are solved numerically to determine the 
nominal controls *^-iQ(e) and 

transcendental nature, the solutions obtained are not 
unique hut depend on the initial guess. The pair of 
nominal controls need to he computed only for the range 
0 5 ® ^ 3-s they repeat after every half orhit. 

(iii) The next step is to choose weighting matrices Q(6) and 
¥(0). In the present analysis equal weights are assigned 
to attitude errors, their derivatives and the control. 

(iv) The matrix-Riccati equation (5.21) is solved backward 
with boundary condition K(® f ■*■“)= 0. In practice, 

0^ = 10 orbits was found to he sufficiently large for the 
solution K(6 ) to reach a steady state periodic character 

with period it . 
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(v) Tile feedback gain matrix 

Ja^ PaM (5.22) 

is then computed in accordance with Equation (5.20), 

(vi) The full nonlinear equations of motion (5.11 ) are finally 
simulated in conjunction with the nominal md the feedback 
control laws synthesized via the linearized analysis. The 
total control surface rotations required during simulation 
are given by 

*l(®) « (5i23a) 

«2(9) = «2o(e) - u(e) (5.23b) 

The earth's shadow criterion (Appendix C) is incorpojrated 
during the simulation. 

(vii) Eor different values of ♦ ,• steps (ii) to (vii) are 

repeated to assess the system performance at different 
times of the year* Note that replacing ♦ by ^ ^ 

Simply leads to an interchange of the nominal control 
histories Simultaneously, if u is 

changed to -u, the governing equations (5»li) i^emain the 
same* Hence, the controller behaviour need be studied 
explicitly only over the half-year interval, 0 5 ^ 1 
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5.4 Results and Discussion 

Figures 5.3 and 5.4 show the nature of the periodic time- 
varying gains for the system. They also indicate the effect of 
the solar aspect angle <p spanning an entire half-year period. 

Typical responses of the system to initial position and 
impulsive disturbances are shown in Figure 5.5 for the case of 
(l> =45®. Figure 5.5a represents a situation when the satellite 
symmetry axis is misaligned from the desired inert ially-fixed 
orientation by 10® in both and cig* controller is able 

to restore the desired alignment in about one orbit. The 
associated feedback control u (Figure 5.5a) and the nominal 
controls (Figure 5.5c) are also indicated. Once the transients 
are removed, the nominal control surface rotations continue to 
counter the gravitational torques and the desired fixed 
inertial orientation is precisely maintained. -The response to 
an initial impulsive disturbance, such as might be caused by 
micrometeorite impacts, is shown in Figure 5.5b. 

The response of the system to the same disturbances 
securing at different times of the year, 4> = 90° and ^ = 135°, 
is presented in Figures 5.6 and 5.7. dlote that the overall 
quality of the response remains essentially the same and only 
some local variations in the response character are observed. 
The system behaviour was found to be quite similar when inertial 
stabilization along other orientations was considered 
(Figure 5.8). 
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The results presented so far correspond to those periods 
of the year when the geost ationary orhit does not encounter 
the earth's shadow. The worst case of the shadow occurs when 
fjf =0 or 1 80® , i . e. , the earth-sun line coincides with the 
line of nodes. IPigure 5.9 represents the system response in 
this situation. SuhseQuent to the demping of the initial 
disturbance effect, the periodic loss of solar control leads to 
a steady-state attitude variation from the desired inertially- 
fixed orientation. The error amplitudes are found to be less 
than 5® under the most severe shadow encounter. Preliminary 
results indicate that through suitable changes in the relative 
weights Q and ¥, the shadow induced errors may be reduced 
further. 

A comment concerning the numerical implementation of the 
feedback gains and nominal controls during the integration 
would be significant here. Since both the system gains and 
the nominal controls are periodic functions, Pourier 
coefficients were fitted to them for convenience of repeated 
generation. Fourier series truncated after seven sine and 
cosine terms first used for both. Although the transients 

damped out effectively, a steady state attitude error (£ 5®) 
remained due to the imperfect cancellation of the gravity 
torque. The nominal controls were then improved through 
linear interpolation over steps of 0.5® . 
practically perfect point ing. ; - 
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The present control policy is valid for a particular time 
of the year. Therefore, continuous updating of the nominal as 
well as the feedback controls would be required as the sun 
position 4 . changes. This would involve a significant amount of 
computation. Considerable simplifications in the software 
implementation would result if the controls developed for a 
particular value of 4> could be used over a range of . This 
would represent a suboptimal control policy. lumerous 
simulations performed in this light showed that the nominal 
controls are quite sensitive to <P and require a more or less 
continuous generation. On the other hand, the feedback gains 
remain effective over a large range of <j> . figures 5.10, 5.11 
and 5.12 show the system performance for <|) = 30°, 60°, 90° where 
the gain histories computed for <J> = 45° have been used. It is 

apparent that a particular gain history remains effective over 
a very large period of the year. Hence, a schedule involving 
only a few gain computations in a year would be necessary. It 
is interesting to compare the suboptimal response for <j> = 90° 
(figure 5.12) with the optimal response for the same case 
presented earlier (figure 5.6). As anticipated, the suboptimal 
case results in a slower rate of damping. The performance 
deterioration, however, is not excessive. 
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5.5 Concluding Reraark: 


The conclusions based on th 
summarized as follovjs s 


e investigation may he 


(i) The potential of a simple two-surface solar controller 
for orienting the symmetry axis of a nonspinning 
spacecraft along any inertially - fixed direction in 
presence of gravitational torques is investigated. 

(ii) A control law consisting of nominal and feedback control 
surface rot-tions has been developed. The nominal controls 
counter the continuous gravity -gradient disturbances and 
the fewclbeck controls impart asymptotic stability to the 
system. 

(iii) Accuratt. attitude pointing requires that the nominal 
controls be updated more or less continuously with the 
time of the year. Hov;ever, the feedback gain scheduling 
need be done only a few times in a year. 

(iv) The control system continues to have a stable operation 
even when tht.* spacecraft passes through the earth’s shadow. 
In shadowu'l orbits, small pointing errors occur which 
depend on the duration of the shadow passage. 



6 , 
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LARGE-MGLE M-TITUDE MAI^Wers np 

JViiKb OP SPACECRAPT SIMMETRI Alls 


The analyses presented so far hare oonsldered the solar 

controller dosip, from the viewpoint of maintaining the 

Spacecraft .ittitude either in earth-pointinff nr ai • 

yuxnxing or along inert ially. 

fixed orientations. It would be of interest to examine the 
feasibility of carrying out arbitrary large-angle attitude 
maneuvers of the spacecraft by means of the two-surfaoe solar 
controller. Such a capability would add versatility to the 
system and rickc it possible to fulfil diverse mission 
requirements. This chapter presents the development of control 
laws for accomplish in,:- large attitude changes of the symmetry 

axis of a nonspinniug spacecraft in the presence of gravitational 
torques. 


6.1 Problem Pomulation and Control Synthesis 

fhe same spacecrait and controller conf igurations as 
presented in Chapter 5 for inertial stabilization are considered 
here. The izovv-vninr, equations of attitude motion, for convenience, 
are restated as 

a!,' cos ag . sin ag - 5(1-I)r^r^ = Q (6.1a) 

°2 + 5(1-I)r r 


= 0 h.j( 6 .j, 62 ) (6.1h) 
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where all syr.bois have the ea.ae meaning ae before. To reeapltu 
late, these e^netlcns present an exceedingly comrlioated, a„n ' 

linear, nonautoncmous system of differential e^atlons. ftrther 

more, they do not possess an eouilihri,™ 

ewrlibrium position and therefore 

appropriate nominal controls are required to n>.+ • 

required to oltain any fixed- 

inertial ori ent at i on . 


The objective of the controller iq +o ^ • n. 

o-Lier IS to drive the symmetry 

axis from an Initial attitude (a., a' a «' \ 

'‘10’ "10’ 20’ “20^ the 

desired final attitiuie(a a„ r,' ) t+ 

^f’ If’ “2f’ “2f^- would, of course 

be deslroblo to achieve the transition in the minimum possible 
time. In .o ncfol, the optimal control would require Independen 
tlme-hlstorus for d^(e) and S^e). However, as pointed out 
earlier, system mechaiisation would be simpler by considering 
equal and opninlto rotations, i.e.,5^(e) = - 62(0). 


for convenience, the governing equations may he rewritten 
in the state vnriaoj ;■ notation. Letting x.^ = a.j, 

X3 = -ind u 6.^ = - 6 ^, they take the form 


X' (6) .1 f(:.:(e), u(0), e) 


( 6 . 2 ) 
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where 


x(0) 

- - r V , V v* T ^ 

1 -.j -2 Xj x^j 

(6.3a) 

f (x,u,e ) 

= ' -^1 ^’2 f, 

(6.3b) 


and the coLrGonc.’rtr- of f are ; 

^1 " '*^2 (6.4a) 

+ 0(1 o.i;- I j onr;; ^ ^ + j cos C 2 1 cos C 2 ) sin u} sec 

(6.4b) 

^3 " '•‘‘4 (6.4c) 


= -3(1-1) 


The 

+ , V 1 

'in Mj] t .i ’ 

*> ,■ 1 cOo - J cos 5 2 1 °cs ? 2 ^ ®cs 

i it s COS 5-| and cos ^2 

u (6.4d) 

given by 

^x ■"= 

(sin X 

1 COB :c-^) CCS 0 - (sin x^) sin 6 

(6.5a) 

= 

(co:"! X 

^ ) cc:; 6 

(6.5b) 


= (sin sin x^) cos e + (cos x^) sine 


(6.5c) 
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00. .in ocs u 

COS ^2='- ^;v ®S ^ ( 6 . 6 b) 

where the co;;.i'rr. cut cf e have the form 

= - r.i:i 4* (sin i cos cos X 3 + cos i sin x^) 

+ ccs<^ (sin cos x^) ( 6 , 7 a) 

.= sin (sin i sin x^) + cos 4 ) (cos x^) ( 6 . 7 b) 

-= “ sin (li (fin i cos x^ sin x^ - cos i cos x^) 

f Cf:;- (P (nin x^ sin Xj) (6.7c) 


i’ht,; ho foimnilated as that of finding u ( 6 ) 

which trihcr. the r.yr.te':. from the given initial state to the 
desired final rtrit. e in minimum time. The resulting two point 
boundary value rrouleni may be solved by techniques such as 
quasilinean cn or variation of extremals. These, however, 
become imi.ract iea.l in l.'nu present situation as the optimality 
condition, 

d n/Bu ^ 0 (6*3) 

where H in the f,;/:rtem i Irani it onian, cannot be solved 
analytical J.y for thci control u in terms of the state and 
costate variables. 

Alternatively, gradient techniques may be used which 

typass the problem of solving the optimality condition (6.8) 
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for the control n. However, the approach requires a knowledge 
of boundary conditions on the costate variables at the 
terminal ti’ e. One of achieving this would be to add a 
penalty function representing a measure of the terminal state 
error to the cost function and letting the final state be 

xf. 

free . Another approach is to require the satisfaction of 
suitably chosen terminal constraints on the states. The latter 
method has boon presently selected to obtain the optimal 
control. 

The control problem is therefore stated as ; 

Find the control u(6) which minimizes the cost function 
®f 

j = J (1) de 

0 

for the system 

x’(6) = f(x(0), u(0), 8) 
with the initial state 

= ‘“10 “io “20 “ 20 '’' 

and the terminal constraints 

T 

where 

4*^ = (1/2) (x^(9f) - “if^ 

= (1/2) (xgCQf) - 


(6.9) 

( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13a) 

(6.13b), 
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t(, j = (1/2) - «2f^^ 

= (1/2) (x^Ce^) -agf)^ 


(6.13c) 

(6.13a) 


6.2 Gradient Technique 

The aesired control may be determined by the approach 
developed by Bryson and Denham^T. the method 

aa adapted to the present problem is briefly described below 


Consider the cost functional 


0 


J =»n(e^), V ^ l.a(e). u(6), 


e ) de 


(6.14) 


for the system 

x'(9) =f(x(e). u(e), e> 

1 • +>, e functional J in the usual manner to 

After augment; mg the i \ a.,, v, 

... equation constraints ( 6 . 15 ),-btoug 
include the differential equ 

Hers 'Me), the differential change 
the use of lagrange multxp variation «u(e). 

dJ due to the differential change j 

may be expressed as 


dJ = fL + + (34>/3® 


® f 

+S + 

0 


X^(3f/3u)l 5 u d0 


(6.16) 


■where X (0) is given t>y 
= - (3Ii/3x) 


— • !E 

^ _ {3 ^ 


(6.17a) 
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with the bomdaay conditions 


X(0^) = (34»/3x) 


(6.l7h) 


Now consider the functional 

ef 


=J (1) 

0 


de 


( 6 . 18 ) 


for the system (6.15). Putting <^ = 0, I = 1 and letting 
p(8) represent the Lagrange multipliers, the differential dJ^ 
is obtained from Equations (6,16) and (6.1?) as 


dJ 


I 

[1 + P^?^J g + J i^(3f/3u) du d6 (6.19) 
^ 0 . 


with 


p’ = - t 3f/3x 1 ^ p 


(6.20a) 


p(0^) =0 

Obviously, Equation (6.20) has the solution p(e) 
Equation (6.19) simply becomes 


dJ = d9^ 
0 1 


Next consider the functionals 


(6.20h) 
0. Hence, 

( 6 . 21 ) 


= «l'i(He^)); i = t,2,5,4 

for the same system ( 6 . 15 )* With ♦ = L = 0 and 


( 6 . 22 ) 
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representing the multiplier vector, one obtains' from 
.Equations (6.16) and (6. 1?) , 


u -1 

■ ' ^ . 0 


•5 


T 


(8f/3u)5ud6 ( 6 . 23 ) 


■with 


= -C 9 f/9x]- 


Ri(ej)= 


(6.24a) 


(6.24b) 


f • 


Eote that dJ. simply equals d4».. Consequently , a change d p. 

-X. -Xp 

in the terminal constraint depends on de^^ and 6u according to 


e. 


dt|>.. = [(9t|;./9x)f] de. + , 

1 10^ f j 1 


J R^(3f/9u) 6u de (6.25) 


Equations ( 6 . 24 ) and (6.25) for i = 1 to 4, may be -written 
in matrix notation as 


R' = - [3f/3x]^ R 


R(e.)= [aV3x] 


T 


(6.26a) 

I 

(6,26b) 


i 

d? =■ { t 3 "^/ax 1 f } p de„ + R^(3f/3u)5u de (6.27) 

®f ^ Jq 


and 



■where 


R = t i R^ ; Rj i R4I (6.28a) 

^ = ('I ' ^2 ^3 ^ 4 )^ ( 6 . 28 h) 

It is apparent that any nominal choice of ^i(9) 

for the system ( 6 . 15 ) 'wotild result in a terminal error 
ir(x(6^)). Corrections d0^ and 6u(0) are now sought which 
change 1jr'(x(0^)) in the direction of satisfying the constraints 
'{»(x(0^)) = 0 and reduce dJ^. One way of doing this would be 
to minimize dJ^ subject to the constraints ( 6 . 27 ) where 
represents the desired change in the terminal error (say 
d^ - - e ilf (x(0^)), 0 < e < 1). Note that dJ^ as well as 

the constraint equation a'it linear in d0^ and 5u(0). To 

ensure that their optimal choices are snail enough for linear 
perturbation theory to be valid, quadratic penalty functions 
in dS^ and 6u(6) are added to dJ to obtain 

0f 

dJ^ = dJ + (1/2)b (de„)^ + (1/2) J ¥(0) (6u)^ de ( 6 . 29 ) 

O O I , Q 

where b is an arbitrary positive constant and ¥(0 ) is a 
positive definite quantity. 


The augmented functional then becomes 

Q ^ 


dJ = de^ + (1/2)b (d0^)^ + (1/2) j"" w(6) (5u)^ d0 

0f.^ 

T r — — ( — 

+ v |{C di|»/3x} f} de^ +J R f/0u)6 


u d0 ■ - di 


?)![ 
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■where v represents the Lagrange multiplier vector. Neglecting 
the change in the coefficients, the differential of dJ be 
written as 

a(aj) = '{ 1 f + b ae .) d(d6.) 

6f 

t (3f/3u)^ Ev+ ¥ 6u] d(6u) dS (6.31) 

Hence, a minimum occurs when : 

d0„ =- (1/b) {1 + [3?/9x } f 1, (6.52) 

1 Of 

6u(8) = - (1/¥) {( 3f/3u)^ Rv} (6.55) 

Substituting Equations (6.52) and (6.33) in Equation ( 6 . 27 ), the 
undetermined constraints v may be obtained from 

dTj> = - (1/b) { (dif>/d0) { (dij»/d0 )'^ V + }g - 

(6.34) 

where „ 

®f 

= J t R^(3f/3u) W"”* (Sf/au)”^ R] de (6.35) 

The solution for v is found to be 

- T 

^ (d4./de) (d4>/d6) X 

■{ df + (1/b) (di/de)}- 

Of 


(6.36) 
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The expected change In is then oWalned hy substituting 
Equ tion (6.32) in liquation (6.21) as 

dJ^ = - (1/b) {1 + [3^ /3xl f } 

When sufficient number of corrections de^ and 6u(e) have 
been applied iteratively, one would expect that df 0 and 
dJ^ 0, which represent the optimal solution. It may, of 

course, be easily shown that the corrections d 0 ^ and 5u(0 ) 
tend to zero as the optimal solution is approached. These, in 
turn, m^ be interpreted as the satisfaction of the transversality 
condition (H|g^ = 0) and the optimality condition (sH/Su = 0) 
associated with the original time-optimal control problem with 
given initial and final states. 

6.3 -Computational Procedure 

The computational procedure may be summarized in the 
following steps : 

(l) Estimate a teminal time 9^ = (8^)^, control history 
u( 0 ) = Uq(9). 

(ii) Integrate the system equations (6.10) with the specified 
initial conditions (6.11). Record i(8), |■(x(9^)) and 

(diT/de ) „ . 

®f 

(111) Integrate Equations (6.26a) backward uith the boundary 
conditions (6.26b), compute the integral I from 
Equation (6.35), and store the history (Sf/Su)"^ R. 
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(iv) Select desired d’J» = - eiji (x(0^)), 0 < ^ < ■]. 

(v) Calculate ^ from Equation (6.36). 

(vi) Compute d9^ and6u(9) from Equations (6.32) and (6.33). 

(vii) Stop if F(x(0^)) = 0 and d0^ = 0 are satisfied to a 
desired accuracy. Otherwise, update 9^ and u(0) 
according to 


new 



+ d0^ 


(6.3Sa) 


u(9) 


new 


= u(9) 


old 


and go hack to step (ii). 


(6.38b) 


6.4 Results and Discussion 

A computer program incorporating the gradient method v;as 
developed. Some computational simplifications were found 
possible without seriously affecting the convergence process. 
The integral ^ was evaluated by assuming the integrand to 
be piecewise constant. The updated 0^ values were rounded off 
to the nearest integral multiple of the step size of the 
integration (0.5°). This makes the storage, updating and 
retrieval of u(0) easily programable. 


The weightages ¥ and b, and the parameter e may be 
varied during the iterative process to improve the 
convergence. Presently, however, constant values ¥ = b = 500 
and e = 0.5 were taken which gave convergence to the desired 
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accuracy (ip^ < 0.005, i = 1 to 4 and } d 6 ^ j < 0.005) in about 
15 - 2o iterations in most cases. 

The optimal control problem was solved for various attitude 
maneuvers over a range of parameters. Typical results 
representing the case of 4» = 45°, i = 23 . 5 °, I = 1/3 and 0 = 10 

are presented and discussed below. Throughout, the initial and 
desired terminal attitude rates are taken to be aero (ct^Q = a^o " ® 
and 

figure 6.1 shows the convergence process for the maneuver 
= - 10°, a 2 Q = - 10° to a-]^ = 25°, a 2 f = “2 

space. As is characteristic of the gradient method, the rate of 
convergence is fast in the beginning and slows down considerably 
as the optimum is approached. The figure also shov;s the initial 
guess and the optimal time history of the control u(6). The 
initial guess of ramp form for u(9) physically represents equal 
and opposite rotations of the control surfaces at a constant rate. 


It appears interesting to examine the effect of different 
initial guesses for ) on the optimal solution, figure 6.2 
shows the results obtained with different ramps selected as u^(6 ) 
with the same guess (9^)^ = 90°. It turns out that different 
optimal controls u^^^(8) are obtained for each case. The 
corresponding trajectories in the , “2 space and the control 
times sre also distinct. This may be interpreted as 

convergence to local stationary solutions in the neighbourhood 




20 40 60 30 1Du 

' , O' 'o , 

Typical convergence to the optimal 
solution for the maneuver iXio=~*^0^ 
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39 “ 

91“ 


0-3 rad /0-5“ 0 
0-2 rad / 0-5" 0 
0*1 rad/O-s'e 


)oDt=9^ 


= 90' 


(0f )opt= 89“ 


(0f )opt = 85 


Fig.6’2 Optimal trajectories and control 

histories obtained with different ra 
type initial guesses for the contro 
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of the initial guesses. Another result with an initial guess 
of rather different character, u^(9) = 0 , is shown in 
Pigure 6.3, 


The influence of the guess (df^Q on the solution is 
indicated in Figure 6.4 for the case of a different maneuver. 
With the control guess ‘WqCq) held the same, the convergence 
to optimal solutions with neighbourhood of 

( 0^)0 is observed. 


The above results clearly indicate the convergence of 
the solution to local optima near the initial guesses ■ 0-^(6 ) 
and Also, it appears that the problem is characterized 

by numerous locally optimal solutions. This may be attributed 
to the highly nonlinear, transcendental nature of the system 
eguations. Nevertheless, the approach is quite useful as it 
enables one to determine a control u (6 ) v/hich can accomplish 
the desired attitude maneuver in inertial space. The transfer 
times obtained are also reasonably small and remain within a 
moderate range despite different starting conditions 
and 

The open-loop control obtained for a specified maneuver 
would drive the spacecraft symmetry axis very close to its 
desired inertial orientation. The satellite attitude, however, 
would deviate from this condition due to the residual terminal 
error and the gravitational torques. Attitude stabilization 










Fig. 6-4 Optimal trajectories and control 
histories obtained with different 
initial guesses for the final time 
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in the desired orientation wuld then he necessary. This may 
he achieved hy applying the nominal and feedback control 
strategies developed in Chapter 5. As an example figure 6.5 
shov/s the attitude transition from = -10°, a2Q = -10° to 
= 250, = 30°. The optimal open-loop control "was 

obtained -with = 0,2 rad/0.5° 6 and = 90° 

(J'igure 6.2). A smooth transition from the given initial 
orientation to the neighbourhood of the desired final 
pointing is achieved by the open-loop control in S9° of the 
orbit. The nominal and feedback controls then stabilize the 
spacecraft along the desired in erti ally -fixed orientation. 
The individual control surface rotations during the entire 
maneuver are : 


for 0 1 ® < 89° 



(6.39a) 


62(9) = 


(6.39b) 

for e > S9° 

61(9) = ^-jQCe) + 

■Uf^(9) 

(6.40a) 


62(9) = ^20^®^ " 

Uf^(e) 

(6.40b) 


■where here denotes the feedback control u(6) of 

Chapter 5. 

The case of yet another arbitrary maneuver (a^Q = -10°, 

= -'^0° to = 15°, “2f = involving the use of 

optimal large-angle transfer and inert ially -fixed stabilization 




Feedback control 


Fig. 6-5 System response and control required 
for large-angle maneuver followred by 
inertially-f ixed attitude stabilization 
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strategies is indicated in figure 6.6. The u 
computed with = 0.1 rad/0.5® e and (8^.) 

(Pigure 6.4). 


(6 ) was 


opt 
„ = 50® 


6.5 Concluding Remarks 

(i) The feasibility of employing a simple two-surface solar 
controller for achieving arbitrary large-angle attitude 
maneuvers of the symmetiy axis of a nonspinning 
spacecraft has been established. 

(ii) An open-loop control strategy for the differential 
control surface rotation is obtained employing the 
gradient procedure. The control law minimizes the time 
required for a specified attitude transition. 

(iii) The numerical procedure was found to converge to 
different solutions depending on the initial guesses for 
the control and the terminal time. With a moderately 
sized controller, the optimal solutions lead to control 
times of a fraction of an orbit for any desired maneuver. 

(iv) After completion of the large-angle maneuver the same 
controller may be used, in the nominal and feedback 
control mode, to stabilize the spacecraft along the 
final ineirtially -fixed orientation. 

(v) The proposed controller adds versatility to the system 
as it makes it possible to achieve different inertially- 
fixed attitudes. 
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7. CliOSIFG OOMMENTS 

7.1 Summary of the Conclusions 

The conclusions based on the preceding investigation may 

be summarized as follows : 

(i) The feasibility of employing a two-surface solar controller, 
representing the minimum hardware configuration, for earth- 
pointing as well as inert ially-fixed attitude stabilization 
has been established. 

(ii) feedback control laws directly governing the control surface 
rotations are synthesized. The resulting control systems 
promise asymptotically stable operation throughout the year. 

(iii) The implementation of the controller for pitch control in 
both earth-pointing and inert ially-fixed modes of 
stabilization is particularly simple, for mult idimensional 
attitude control, implementation of time-varying gains is 
involved. Considerable computational simplifications may 
be achieved through the use of suboptimal control policies. 

(iv) The control systems, using moderately sized controllers, 
are found to be capable of stabilizing the spacecraft 
even along gravitationally unstable orientations. Time 
constants of the order of a few orbital degrees ar^ 
possible and may be reduced further through the use of 
increased gains. 
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(v) The effectiveness of the controller in maint aining 
inertially-fixed attitudes for a nonspinning spacecraft 
in the gravity -gradient field suggests new applications 
for the concept of solar pressure control. 

(vi) An open-loop control law for achieving arbitrary large- 
angle attitude maneuvers is established. This capability 
makes the system more versatile and opens the possibility 
of designing satellites for diverse missions. 

(vii) Due to the loss of solar control during encounters v/ith 
the earth's shadow, small attitude deviations are 
observed. On emergence from the shadow, however, the 
controller is able to reorient the satellite precisely 
along the desired orientation, 

(viii) An analysis accounting for the possibility of one of the 
control surfaces being shadowed by the spacecraft would 
be highly dependent on the relative shapes and sizes of 
the satellite and the controller. Preliminaiy results 
suggest the effect to be of little consequence for earth- 
pointing applications and during inertial stabilization 
along orientations that are not too close to the sun-line, 
(ix) Attitude control systems have been presented for spacecraft 
in circular orbits. Through minor modifications of the 
control policies, solar controllers applicable to 
satellites in elliptic orbits may be developed. 
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(x) The concept of solar presstire control does not involve 
any mass expulsion schemes and/or active elements 
requiring large power consumption. The semipassive 
character of the control systems promises increased life- 
times for the spacecraft. 

7.2 Recommendations for Puture Work 

On the basis of the present study, the following 

recommendations are made for future work : 

(i) A ^stematic study may be carried out to determine the 
best orientation of the rotation axis of the two-surface 
controller relative to the spacecraft body. The analysis 
would likely involve optimizing the control torque 
components averaged over the year as a function of the 
orbital inclination. 

(ii) The inertial orientation problem may be extended to three- 
axis control through the use of three rotatable control 
surfaces. 

'iii) The concept of nominal control may be applied to counter 
the continuous disturbance due to radiation pressure 
arising for spacecraft with a nonzero offset between the 
centers of mass and pressure. 

(iv) The possibility of modifying the approach for solar control 
of spacecraft having a limited attitude determination 
capabilities could be explored. 
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(■v) Utilizatation of solar radiation pressure for attitude 
control has so far been limited to rigid spacecraft. On 
the other hand, increasingly flexible configurations are 
being contemplated to meet the challenges of the future, 
feasibility of the concept of solar control for flexible 
spacecraft represents an exciting area for future 
investigations. 

(vi) The development of solar controllers for achieving 

controlled variations of the orbital elements and orbital 
transfers represents another field with important 
challenging problems. 
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iPPENDK A 

The components of the earth-sun unit vector e in the 
satellite principal axis system are given hy : 

= sin [ - sin i cos y cos 6 + cos i(sin 6 sin Y cos ^ 
- cos 0 sin 6)1 

+ cos <|> [ sin 0 sin 0+ cos 6 sin y cos 0} (A-1 ) 

e = sin <|)[ sin i(sin y cos X - cos y sin B sin X ) 

if 

+ cos i { sin 9 (cos Y cos X + sin Y sin 6 sin X ) 
+ cos 9 (cos B sin X ) }] 

+ ' cos (|i I - sin 9 (cos B sin X) 

+ cos 9(cos YGos X + sinY sin B sin X) 1 (A-2) 

= sin t - sin i(sin y sin X + cos Y sin B cos X) 

+ cos i { sin 0 (sin y sin boos X - cosy sin X ) 

+ cos e (cos B cos X )}1 
+ cos #[- sin 9 (cos B cos X) 

+ cos 6 (sin Y sin B cos X - cos y sin X )1 (A-3) 
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The functions , f ^ appearing in Equations (3.11) 
are given by : 

^ ~ ^ ** ^ ~ ^ ^ ^ ^ ^ 

-Cly Gos^ ^2 sin^ X)cos y 

-Cly - Ij.) sin Y sin 3 sin 2 X 1 g’ 

'<• { (1/2) (ly - l2)(2 sinY cos geos 2X 

- cos Y sin 2 8 sin 2X ) I siny cos 3 1 X* 

Jv 

+ (1/2)(I^ - I^) 8’ ^ sin 3 sin 2 X 

“ ^ “ ^2 ) y'P' sin 2 8 

4 1 Ijj. - (ly - I^)cos 2 X] 3’ X' cos 8 

- (ly - ^2^ 

- (1/2) [I sin 2 y cos^ 8 + { (L sin^X + I cos^X)sin^ 3 

A y z 

- (ly cos^X + sin^X )} sin 2 y 
-)■ (I - I ) cos 2Y sin ^ sin 2 X 1 

tj ■“ 

“ (3/4) I (2I - I - I_) sin 2 y cos^3 

A y 2 

2 

- (ly-I^) (sin 2 y cos 2X + sin 2 y sin 3 cos 2 X 

- cos 2Y sin 3 sin 2X )) 


(A-4) 
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fs - - 008 28 + J. Sloe ) 008 Y 

+ V^,) (y 8lnB-y) 3^ 2 ;^ 

+ t(I, - ^ 


ly sin X + oogSj J 
+ (ly 008^ 1+ I_ 31^2 

- 


^ )] Y 


COSY 


A QOS Y sitt 6 cos 2 X 


+ (1/2) (I 


x~Iy sin X - cos X )(y' 2 _ CQg2Y 


+ 3 sin^ Y ) sin 2 g 

+ (ly - I 2 )(B' X' _ sij:i 2Y cos 6) sin 2 X 
{ + (l^ - I^) cos 2X] yl X' cos B 


(A-5) 


^X = y Y ' (sin Y + e ' ) cos 6 + B r cos Y sin B } 

- (V2) (I^ - l^) {(g. ^ gi^ 

- (y' cos B + cos Y sin B)^ 

+ 3(sin^Y sin^ B _ cos^ Y )] sin 2 X 
(5/2) (ly. - l^) sin 2Y sin 0 cos 2 X 


(A-6) 
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APPEH)IX b 

The quantities (i = 1 to 5) appearing in 

Equations (5.15) are given hy . : 

= (1 - I)k/cos « 2 e (B-l ) 

a2 = -3(1 - I) l^/cosa2g 

- (20/ cos 02 ^) t I q^l r^ sin 6 ^q - | q 2 | ^2 sin 6 2^1 

(B-2) 

= -3(1-1) I 2 /C 0 S 

-(2G/gos a2e^ ^ k-ll sin - I q2i S2 sin 5 go I 

(B-5) 

a^ = -(2C/cos U-jltt^ sin + (q^/2) cos 5 ^q] 

(B-4) 


a^ = (2C/C0S I q2l r'^2 *^20 ^20 ^ 


= -(1-I)k cos 

• 

(3-5) 

(3-6) 

Eg = 3(1-I)m^ 

-2GClq^l3:^ cos - 

1 rg cos 'SgQl 

(B-7) 

= 3(1-l)m2 

-20[ Iq^! cos - 

1 qgl Sg cos 5 20 1 

(B-S) 

= -2G j q^l [ t^ cos 

- (q^/2) sin «-,q1 

(B-9) 
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^5 " ft2 cos 620 - (q2/2) sin (B-lo) 

where 


11 " ^“le °°s“2e^ + (1/2) (sin sin 

(B-11) 

1 2 = -(1/2)1 (sin 2 a^^ giK a,^^)cos^9 + (cosa^^ cosa 2 ^)sin 20] 

(B-12) 



= ( 1 / 2 ) [ (sin 

sin 

2 “ 2 e) 

2 

cos 

6 





+ ( COSO ^ ^ 

cos 
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cos ®-j 

^ sin ' 
e 

“ 2 e 
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le 
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^30 
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COS 
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Figure 0-1 sho-ws the geometry of the earth's shadow which 
is assumed to he an infinite half-cylinder. The angle between 
the earth— sun line and the local vertical is given by 

cos = (1/R)R . e (0-1) 

The vectors R and e may be written in terms of their components 
along X,Y,2 axes as 


R = (R COS'S ) J + (R sin 6) K (0-2) 

e = (-sin<I> sin i) T + (cos<|) ) J + (sin<J) cos i) f 

(0-3) 

Substitution of Equations (0-2) and (0-3) in Equation (0-1) 
yields 

cos = cos 4> cos 0+ sin(j> cos i sin 8 (0-4) 

The distance of the satellite from the earth-sun line is given 

R = R I (1-cos^’J' )^ 1 (0-5) 

s 

Ihe condition for the satellite to be in the earth's shadow 
is then readily expressed as 






0 


cos f < 0 


and 


(0-6) 






Fig. C-1 Geometry of the earth's shadow 
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It may te easily shovm that the satellite is eclipsed 
by the earth for the longest period whai $ = 0 or 180°. 
The longest shado-w passage for the case of a geostationaiy 
orbit occurs for about 18° of the orbit. 
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